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Abstract 



The discovery of supersymmetry is one of the major goals of the current exper- 
iments at the Tevatron and in proposed experiments at the Large Hadron CoUider 
(LHC). However when sparticles are produced the signatures of their production will 
to a significant degree depend on their hierarchical mass patterns. Here we investigate 
hierarchical mass patterns for the four lightest sparticles within one of the leading 
candidate theories - the SUGRA model. Specifically we analyze the hierarchies for 
the four lightest sparticles for the mSUGRA as well as for a general class of super- 
gravity unified models including nonuniversalities in the soft breaking sector. It is 
shown that out of nearly 10^ possibilities of sparticle mass hierarchies, only a small 
number survives the rigorous constraints of radiative electroweak symmetry break- 
ing, relic density and other experimental constraints. The signature space of these 
mass patterns at the LHC is investigated using a large set of final states including 
multi-leptonic states, hadronically decaying rs, tagged b jets and other hadronic jets. 
In all, we analyze more than 40 such lepton plus jet and missing energy signatures 
along with several kinematical signatures such as missing transverse momentum, ef- 
fective mass, and invariant mass distributions of final state observables. It is shown 
that a composite analysis can produce significant discrimination among sparticle mass 
patterns allowing for a possible identification of the source of soft breaking. While 
the analysis given is for supergravity models, the techniques used in the analysis are 
applicable to wide class of models including string and brane models. 
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Chapter 1 
Introduction 



While the Standard Model of particle interactions is highly successful, important 
gaps remain in extending the model to a more complete unification, including the 
electroweak and the strong, and eventually the gravitational interactions. Over the 
past decades, supersymmetry (SUSY) has turned out to be one of the leading can- 
didates for physics beyond the Standard Model. In this thesis, we investigate the 
signatures at the CERN Large Hadron Collider (LHC) for some of the supersymmet- 
ric theories. This analysis can be thought as a map from the parameter space of the 
underlying theories onto the signature space of the LHC, as indicated in Fig. (11. ip . 
The parameter space is spanned by the input parameters of the theoretical models, 
and the number of the independent input parameter can be as many as ~ 110 in the 
Minimal Supersymmetric Standard Model (MSSM). It is not realistic to investigate 
models with such a large parameter space for their LHC signatures. However, it is 
possible to do so in some well motivated models where the dimensionality of the pa- 
rameter space reduces significantly. Specifically we will focus on models where the 
dimension M of the parameter space is small, as often M = 4. 

The strategy of this analysis is as follows. We start with a point in the M- 
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Figure 1.1: The map between the parameter space and the signature space. 

dimensional parameter space, investigate the theoretical predictions of this model 
point, simulate its behavior within the Large Hadron Collider, and finally extract 
many possible LHC signatures which is represented by a point in the N-dimensional 
signature space as shown in Fig. (11. ip . The dimensionality of the signature space 
can be expanded by adding more LHC signatures. However, many of the signature 
channels are strongly correlated so the independent set of signatures is typically much 
smaller than the number of signatures investigated. 

The analysis in this thesis focuses on supergravity models, including the min- 
imal supergravity grand unification models (mSUGRA) and SUGRA models with 
nonuniversalities (NUSUGRA). The mSUGRA model depends on four soft breaking 
parameters, and the SUGRA with nonuniversalities that we investigate here contain 
two more which characterize the nonuniversalities in different sectors. We perform an 
exhaustive scan with Monte Carlo simulations in the mSUGRA soft parameter space 
(4-D parameter space) and in NUSUGRA parameter space (6-D parameter space). 
The model points that pass the various cosmological and collider constraints are clas- 
sified by their mass hierarchical patterns. These hierarchical patterns are defined by 
the four lightest sparticles. We will see that the hierarchical mass patterns to a great 
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degree influence the LHC signatures. Most of the analysis given here will attempt to 
correlate the hierarchical mass patterns with specific signatures at the LHC in the 
mSUGRA model. Extension of the analysis to NUSUGRA will also be discussed. 

The thesis is organized as following. In chapter ([2]) we give a brief introduction 
to the Standard Model (SM), and in chapter ([3]) we discuss Supersymmetry (SUSY), 
and the Minimal Supersymmetric Standard Model (MSSM). In chapter (jlj), we give 
an introduction to supergravity unified models. We list the various experimental 
constraints that are imposed on the supergravity models in chapter ([5]). 

An analysis of hierarchical mass patterns within the mSUGRA and NUSUGRA 
models is given in chapter ([6]). Here the correlations between the sparticle patterns 
and the nature of the soft breaking are also analyzed. In chapter (JTj), we give a detailed 
description of our simulations of the CERN Large Hadron Collider which includes 
the various steps of our simulation, the detector cuts, and the various signatures 
we investigate. In chapter ([8]), we analyze the signatures arising from the sparticle 
patterns, and also use the signatures to discriminate the patterns. Several important 
kinematical signatures are analyzed in chapter ([9]), including the missing transverse 
momenta, effective mass, and the dileptonic invariant mass distributions. In order to 
utilize as many signatures as possible, a global analysis using the so called fuzzy vector 
technique is given in chapter ( ITOi) . An analysis regarding the signature degeneracies 
from different models and the ability of resolving the parameter space using the LHC 
data is carried out in chapter ( fTTl) . 

In addition to the signature analysis for the LHC, we also investigate other exper- 
imental signatures within the context of the sparticle pattern discrimination. These 
include the analysis of Higgs production at Tevatron in chapter ( fT2i) . the Bg fi^fi^ 
constraints in chapter f|T3l) . and the direct detection of dark matter in chapter f|T4|) . 
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Chapter 2 

The Standard Model 



The Standard Model (SM) [1, 2J of particle physics is a theory that describes three of 
the four known interactions, which are the electrodynamics, the weak interactions and 
the strong interactions. The SM is based on the gauge group SU (3)cxS'f/ (2)^ xf/(l)y 
where C stands for color, L stands for left chiral, and Y stands for hypercharge. There 
are three generations of quarks and leptons in the SM which are represented by left- 
handed doublets and right handed singlets in the SU{2)l gauge group. 



Ul 
ydLj 



URi, 



dm', 



em 



(2.1) 



where ^l,r = Pl,r^, Pl = {I - 75)/2, Pr = {I + 75)/2, and ^ = 1,2,3 is the 
generation index. In order to give mass to quarks and leptons, as well as to vector 
bosons, an additional scalar SU{2)l doublet is introduced in the theory 



(2.2) 
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The three fundamental interactions (the electroweak and the strong interactions) are 
mediated by eight SU{S) color gluons, G^; three SU{2)l gauge bosons, Ay, and one 
U{1) hypercharge gauge field, B^. All the above gauge bosons are realized in the 
adjoint representations of their corresponding gauge groups, and the strength of the 
interactions are described by their coupling constants g^, g2 and g'. 

The dynamics of the Standard Model consists of the following three parts: 

1. Gauge interactions, 

2. Yukawa interactions, 

3. Higgs potential. 

The gauge interactions arise via the gauge covariant derivative 



D„ = di,-i 



a=l 1=1 



(2.3) 



where T^=(A"/2; 0) for (quarks; leptons, Higgs), where the A" are the eight Gell- 
Mann matrices; Tl={a^/2; 0) for SU{2)l (doublets; singlets), where cr* are the three 
Pauli matrices; and Y is the hypercharge defined by Q = T| + Y/2. 

Since the explicit fermion mass terms violate the gauge symmetries, adding mass 
terms in the SM is done through the Yukawa interactions and the mechanism of 
spontaneous symmetry breaking. Thus we consider the Higgs potential given by 

y (0) = + A((/)t0)2 (2.4) 

where /x^ and A are positive. The tachyonic mass term for the Higgs field gives rise to a 
spontaneous symmetry breaking which results in a non- vanishing vacuum expectation 
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value (VEV) of the Higgs field 



(0)c 



( ^ 



(2.5) 



where v = fi/^/X. The neutral Higgs gets redefined with respect to this VEV, 
(v + K)l \pl and the new Higgs boson possesses a positive tree level (mass)^ of M| 
2/i2 > 0. 

The Yukawa interactions which preserve the gauge symmetries are given by 



Vy = Xf^Uen, + X^^^UuRj + Xf^q^dn, + h.c. (2.6) 



where = ia'^cf)* is the complex conjugate of the Higgs doublet given in Eq. (12. 2p . 
and Aj-^'"''^^ are the Yukawa coupling constants where i,j = 1,2,3 are the genera- 
tion indices. The Yukawa interactions are responsible for giving rise to the fermion 
masses. The flavor basis of the fermion species are not necessary the same as the mass 
eigen states of the particles, and the unitary transformations between the flavor eigen 
states and the mass eigen states is accomplished via the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. 

The spontaneous symmetry breaking also rotates the four SU{2)l x f/(l)y gauge 
bosons to their mass eigen states by means of the gauge interaction term of Higgs 
fields, {A^,A^} {W^,W^} and {A^,B} {A, Z} where the photon field A re- 
mains massless. Introducing the weak mixing angle, 9w defined by tanOw = g'/g, we 
can express the photon and Z boson fields as follows 

A = cos 9wB + sm9wA\ (2.7) 
Z = - sinOwB + cosewA^. (2.8) 
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The and Z boson masses at tree level are given by 

M„ = i^ = -%M„ M. = ^, (2.9) 
2 2V2A cosOw 

showing that the Higgs mass sets the electroweak mass scale. 

Standard Model is a highly successful model. It predicted the existence of the W 
and Z bosons before these particles were observed. Since then, the SM has passed a 
variety of experimental tests with great accuracy. However, from the theoretical side, 
many aspects of the Standard Model are not satisfactory. There is no real unification 
among the electroweak and the strong interactions as the SM gauge group is a product 
of three different gauge groups. 

However, the most serious problem within the Standard Model seems to be the 
"gauge hierarchy" problem, or "naturalness" problem which concerns the renormaliza- 
tion of the Higgs mass. Within the SM, the one loop corrections to are quadratic 
in the momentum cutoff. If the cutoff momentum were to placed at the GUT scale 
Mg ~ 10^^ GeV, the tree level Higgs mass is forced to the same scale, in order 
to reproduce the electroweak scale. This then results in the "fine tuning" problem 
because of the huge difference between the GUT scale and the electroweak scale. 
Another shortcoming of the Standard Model is that gravity is not explained in its 
framework. 
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Chapter 3 

Supersymmetry &: the Minimal 
Supersymmetric Standard Model 



Supersymmetry (SUSY) [3J is a symmetry between bosons and fermions, i.e. it re- 
quires that the number of bose and fermi hehcity states in a multiplet be equal. 
A supersymmetry transformation turns a bosonic state into a fermionic state, and 
vice versa. The generator Q and its hermitian conjugate of a supersymmetry 
transformation obey the so called "Graded Lie Algebra" which includes the following 
commutation and anti-commutation relations 



{Qa,Ql} = 2<P^, (3.1) 



{Qa.Qp} = {QlQ\} = 0, (3.2 



[Qa,P,] = [QlP,] = 0, (3.3) 
[P,,Pu] = 0, (3.4) 

where cx^ = (—!,—(?) with a' being the Pauli matrices, and the undotted (dotted) 
indices, a = 1,2 (d = 1,2) are introduced when a four- component Dirac spinor is 
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decomposed into two two-component Weyl spinors. 

The simplest SUSY multiplets are the massless chiral multiplet and the vector 
multiplet. The left chiral multiplet consists of one left handed Weyl spinor (spin 
1/2) and its superpartner, one complex scalar field (p (spin 0, left-handed). The Weyl 
spinors can be used to represent fermionic matter fields such as quarks and leptons 
and the scalar fields are the superpartners, i.e. "squarks" and "sleptons". The vector 
multiplet in the Wess-Zumino gauge consists of one vector field (spin 1) and one 
Majorana spinor (spin 1/2). The vector fields represent the gauge bosons, and the 
additional spinors arc the superpartners of the gauge bosons called "gauginos" . 



quark 



Ul 
dL 



URi 

dRi 



squark 



.dL 

URi 

dRi 



lepton 



l^L 



slepton 



l^L 
.^L 



e-Ri 



GRi 



Higgsino 



Hi 



Ho 



Hf 
H{ 



Higgs 



Hi 



H, 



H', 
Hr 
Ht 



Table 3.1: Chiral supermultiplets in the Minimal Supersymmetric Standard Model. 



gluon 


G^,(a = l,...,8) 


gluino 


G", (a = 1,...,8) 


SU (2) gauge boson 


^i.,(^ = l,2,3) 


SU (2) gaugino 


AS (^ = 1,2, 3) 


U{1) gauge boson 


Bl 


U{1) gaugino 





Table 3.2: Vector supermultiplets in the Minimal Supersymmetric Standard Model. 



The Minimal Supersymmetric Standard Model (MSSM) is the simplest supersym- 
metric extension of the Standard Model. One promotes each of the Standard Model 
particles to either a chiral or a vector multiplet, which makes the particle content 
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roughly twice as big as in SM as shown in Table (13. 11 13.21) . In the Higgs sector one 
has two Higgs doublets, one of which (-^2) gives mass to the up quarks and the other 
(Hi) gives mass to the down quarks and the charged leptons. 

In total there are 32 supersymmetric particles. These include 4 Higgs boson states, 
of which three {h, H, A) are neutral, the first two being CP even and the third CP odd, 
and one charged Higgs H"^. In the gaugino-Higgsino sector there are two charged mass 
eigenstates (charginos) Xili,2) ^^ur charge neutral states (neutralinos) Xi'=i,2,3,4; 
the gluino g. In the sfermion sector, before diagonalization, there are 9 scalar leptons 
(sleptons) which are superpartners of the leptons with left and right chirality and are 
denoted as: {cl^r, P'L,r, t'l,r^ ^sl^ ^iml^ ^tl)- Finally there are 12 squarks which are the 
superpartners of the quarks and are represented by: {ul,Ri cl,r, iL,R, dL,R, sl,r, 
Mass diagonal slepton and squark states will in general be mixtures of L, R states. 

In MSSM the superpotential with R-parity conservation is given by 

W = U^'Y^QHu + D^'YaQHci + E^vMd + I^H^H^ (3.5) 

where Yu^d^g are matrices in family space. One can add a large number of soft terms 
to the Lagrangian. Examples of R-parity conserving terms are 

Csoft = Cfoft + '^io/r (3-6) 
The soft SUSY-breaking Lagrangian contains scalar couplings 

= u^KQH^ + d^'hdQHd + -e^'KLHd + BH^Hd + h.c. (3.7) 
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where hu^d,e are 3x3 matrices. There are also scalar masses 

+vF^ml vF + dP^ml dP + e^W? (3.8) 

where M| , M|, m|, m|, and m| are 3x3 matrices in family space. More generally, as 
already noted, there can be as many as ~ 110 parameters in the soft terms. However, 
as we will see, the sets of parameters decrease significantly in Supergravity Unified 
Models. 
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Chapter 4 

Supergravity Unified Models 



Supersymmetry provides a solution to the gauge hierarchy problem, which makes it 
an attractive candidate for the new physics beyond the Standard Model. The main 
hurdle in the development of realistic supersymmetric models in the early days was 
the difficulty of breaking supersymmetry in a phenomenologically viable manner. In 
the framework of supergravity, this problem is solved by the inclusion of the gravity 
into the analysis, which promotes supersymmetry from a global symmetry to a local 
symmetry [HE], [6]. 

To construct viable supergravity models one must couple = 1 supergravity 
with vector gauge fields and with matter. This construction is often referred to as 
"applied supergravity" [SI [7l E] and herein one couples = 1 supergravity with 
= 1 Yang-Mills fields in the adjoint representation of the gauge group G (where G 
could be SU{3)c x SU{2)l x U{1)y, or SU{5), SO{10), E{Q) etc), and with N = 1 
matter which contains quarks and leptons and Higgs fields which belong to anomaly 
free combinations of representations of the gauge group. The most general effective 
Lagrangian thus constructed depends on three functions which are the superpotential 
W{(j)i), the Kahler potential K{(j)i,(f)l), and the gauge kinetic function faf3{4'i,<Pl) 
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where a, (3 are adjoint representation gauge indices, 0j are the spin zero components 
of the left handed chiral muhiplet consisting of Xi)- We note that W, K and f^p 
are hermitian. In fact W and K enter in the effective theory only in the following 
fixed combination 

g = K^K + ln[K^WW^] (4.1) 

where 

K = 1/Mpi (4.2) 
and Mpi is the Planck mass defined in terms of the Newton's constant Gjy hj 

Mpi = (SttGtv)"^/^ = 2.4 X lO^^GeV. (4.3) 

The above implies that the supergravity Lagrangian effectively depends only on two 
functions which are fajj and Q. From the above one can easily check that the effective 
theory is invariant under the Kahler transformation as given below 

K^K- /(0,) - /t(0,), W ^ e^'^iy. (4.4) 

It is useful to introduce the so called Kahler metric as follows 



K] = K^^ ^ = K-'g^;. (4.5) 



The case K = - (j)i(j)\ gives i^j = 5] which is referred to as the flat Kahler metric. 

One of the most important results that emerges from the applied supergravity 
analysis is that the effective potential of iV = 1 theory takes the form [HI E] 

V = K-^e-^ [{g'^)]g,g^' - s] + C [R<rXp\ D^Df,. (4.6) 
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In the above g is the gauge couphng constant, and (/ and {Q are the matrix 
inverses of fap and Qj while is given by 

Da = K-^g^\T'')ijZj (4.7) 

where is the group generator. An alternative form which is often useful is to write 
the scalar potential explicitly in term of W and K, and one then has 



kK 



+ Vd (4.8) 



where Vd is as given as before. In the applied supergravity construction the kinetic 
energy of scalar fields is given by 

- K':^{D^ct>,) {DM^ (4.9) 

where is the gauge covariant derivative. 

Before proceeding further we comment on the // term that arises in the Higgs 
bilinear term in the form fiHiH2 in the superpotential. For phenomenological reasons 
II must be of electroweak size, and thus one might speculate on the origin of this term 
in the superpotential. In fact it is not difficult to see how such a term can arise. The 
simplest way to envision the generation of such a term is via the Kahler potential. 
Thus one can write the Kahler potential for the MSSM case so that 

K ^ Ko{(j)^(j)l) + coH^H^ (4.10) 

where 4>i are the MSSM scalar and Cq is dimensionless. Note that 11QH1H2 is the most 
general bilinear term which one write without introducing dimensioned parameters 



23 



in the Kahler potential. Specifically operators of dimension greater than 2 will be 
suppressed by 1/Mpi and thus their contributions will be small. Next one can make 
a Kahler transformation and move the C0H1H2 term from the Kahler potential to the 
superpotential [9]. Indeed a term of this type arises naturally in string constructions 

m- 

After the transformation the superpotential has the form 

lye^'/ = W + coK^WHiH2 + ■■ (4.11) 

As is discussed below spontaneous breaking of supersymmetry by gravity mediation 
gives a non- vanishing VEV for W and we define 

fio = Cot^'{W) (4.12) 

where as will be seen below the quantity /^"^{W) is of electroweak size. 

We now turn to the issue of breaking of supersymmetry. One of the reasons that 
globally supersymmetric models do not have phenomenologically acceptable breaking 
of supersymmetry is that here one has a positive definite potential which after spon- 
taneous breaking of supersymmetry leads to a large non-vanishing vacuum energy. 
This problem is corrected in supergravity models. Here as seen in Eq. (14. Sp one finds 
that the potential contains a term with a negative sign and thus after spontaneous 
breaking the vacuum energy can be fine tuned to zero. 

The central assumption of breaking of supersymmetry in supergravity models is 
that supersymmetry is broken in a hidden sector and the breaking is transmitted by 
the gravitational interactions to the visible sector. A specific illustration of this comes 
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about as follows: One writes the superpotential in the form [6], [TT] 



(4.13) 



where PVvis contains fields of the visible sector which are the MSSM fields including 
quarks and leptons and Higgs fields, and Whid contains fields in the hidden sec- 
tor where supersymmetry breaks. The breaking gives a non-vanishing VEV so that 
(W) = (PVhid)- The size of (PVhid) is estimated to be rn?M^i where m is an interme- 
diate scale so that m ~ 10^^ GeV. The breaking is transmitted to the visible sector 
by gravitational interactions producing soft breaking terms. The gravitino develops 
a mass which is M3/2 = K^e^^'^\W\. Similarly the scale of soft breaking that enters 
the scalar sector is given by 



Thus with m ~ 10^^ one finds mo ~ 10^ GeV, i.e., of electroweak size. Other soft 
terms can be generated in a similar way. A remarkable aspect of the analysis is that 
in supergravity grand unified models the soft breaking is independent of the grand 
unified unification scale Mq which cancels out in the low energy theory [6l [121 [13] 

The gauge kinetic energy terms that generate masses for the gauginos are exhibited 
below 



In general the gauge kinetic energy function /q,^ has a non-trivial field dependence 
involving fields which transform as a singlet or a non-singlet irreducible representation 
of the underlying gauge group. After the spontaneous breaking of supersymmetry the 




(«:'(W^hid))'. 



(4.14) 



■gauge 
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above lead to gaugino masses. If one assumes that the fields transform as singlets 
of the underlying gauge groups, then the gaugino masses at the GUT scale will be 
universal and generate a term of the form mii2\a^a- Splitting of the gauge masses 
can be obtained by the assumption that fap have fields which transform as non-singlet 
irreducible representation of the underlying gauge group. 

The phenomenology of supergravity (SUGRA) models has been discussed since 
the inception of these models (for reviews see [HI [151 IIHl [U] ^-^id there exists now 
a considerable amount of literature regarding the implications of SUGRA (for early 
works see [H], for more recent works see [HI [201 [2l] and [221 [231 [2^1 [251 [26]. for 
works with nonuniversalities see [27], and for works with hierarchical breaking and 
with U{1) gauge extensions see [28l [291 [30]). While many analyses of the mSUGRA 
parameter space have been limited to the case of vanishing trilinear couplings, several 
recent works [3ll [32l [331 [311 [201 [211 [351 [361 [37] have appeared relaxing this assumption, 
and new portions of the parameter space have been found consistent with all known 
experimental constraints on the model. 

As mentioned already in the previous chapter, there are 32 different supersym- 
metric particles, or sparticles. If all the 32 sparticle masses are treated as essentially 
all independent, aside from sum rules (for a pedagogical analysis on sum rules in the 
context of unification and RG analysis see [3H]) on the Higgs, sfermions, chargino 
and neutralino masses, then without imposition of any phenomenological constraints, 
the number of hierarchical patterns for the sparticles could be as many as 0(10^^) 
or larger. This represents a mini landscape in a loose way reminiscent of the string 
landscape (which, however, is much larger with as many as O(10^°°°) possibilities) 
[39]. (Here we refer to the landscape of mass hierarchies and not to the landscape 
of vacua as is the case when one talks of a string landscape. For the string case the 
landscape consists of a countably discrete set, while for the case considered here, since 
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the parameters can vary continuously, the landscape of vacua is indeed much larger. 
However, our focus will be the landscape of mass hierarchies.) Now, the number 
of possibilities can be reduced by very significant amounts in supergravity models 
with the imposition of the constraints of radiative electroweak symmetry breaking 
(REWSB) which we discuss below. 



4.1 Radiative Electroweak Symmetry Breaking 

In the minimal supergravity unification (mSUGRA) the potential at the GUT scale 
which gives rise to the soft breaking of supersymmetry is given by 

F = 5^ffp + m^5^0l0, + Aol^(3)+i?oW^(^) + mv2 Yl (^.15) 

a i a=3,2,l 

where 0j are the scalar fields, Xa{<y = 3,2,1) are the gauginos corresponding to 
SU{3)c, SU{2)l and f/(l)y gauge groups, Aq is the coefficient of the trilinear cou- 
pling, and Bq is the coefficient of the bilinear coupling. For MSSM, VT^^) takes the 
form 

1^(2) = ^.eijHlHi (4.16) 
where i,j = 1, 2, and W^^^ takes the form 

W^'^ = >!i;\.H,u^ + ASf g.i/icif + A^^.i/ief (4.17) 

where A^^J''^'^'' are the Yukawa couplings (analogous to those in Eq. (12. 6p ) and Hi and 
H2 are the two Higgs doublets. 

The soft breaking for the minimal supergravity model is thus given by the param- 
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eters 



mo, mi/2, Ao, Bq. 



(4.18) 



While the mSUGRA model is initialized at the GUT scale, the experiments are car- 
ried out at the electroweak scale. One needs the Renormalization Group Equations 
(RGE) to connect these two domains [lO]. For discussion of the electroweak symme- 
try breaking, one needs the effective Higgs potential which is V = Vq + AVi where Vq 
is the tree part and AVi is the one loop corrections [ITl H2] 



Vo = mi\H,\' + mi\H2Y - mi{H,H2 + h.c.) + -{g^ + 9p)im' - \H2\y, (4.19) 



and 



M2 



e3/2Q2 



(4.20) 



Here, mi{t), g2{t), grit) are all "running" parameters at scale Q where t = ln(M^/(5^). 
Thus US] 

m^{t)=mUt)+i^\t), 1 = 1,2; (4.21) 



with the boundary conditions at the GUT scale Q = Mdt = 0): 



(4.22) 



(4.23) 
(4.24) 



The coupling constants are unified at the GUT scale too 



"2(0) = (5/3)ay(0) = ac. 



(4.25) 
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In Eq. (14.11) . Ma = Ma{vi,V2) is the tree level mass of particle a as functions of the 
Higgs VEVs, Vi = {Hi), and Sa and are the spin and number of helicity states of 
particle a. 

At the GUT scale, all scalar particles are of positive mass^ values ttiq. When one 
integrates down in the energy scale with RGEs, some mass^ values can turn negative 
which will break down the SU (2) x f/(l) symmetry as we have seen in the chapter ([2]). 
The SU{2) X U{1) symmetry breaking requires two conditions: (1) the determinant 
of the mass^ matrix be negative so that there exists one negative eigenvalue; (2) the 
potential be bounded from below. We apply these two conditions to the tree level 
Higgs potential and get 



The above two conditions cannot be satisfied simultaneously at the GUT scale since all 
scalar particles have the same mass. However, when one evolves down from the GUT 
scale to the electroweak scale, the heavy top quark contributions to naturally 
pushes it to turn negative so that electroweak symmetry breaking can be achieved. 
The fact that the top quark must be heavy was one of the predictions of the su- 
pergravity theories I^^Sj (see [44J for a recent review of radiative breaking). Because 
the electroweak symmetry breaking is driven by the quantum loop corrections, this 
mechanism is therefore known as radiative electroweak symmetry breaking (REWSB). 

Turning now to the minimization of the Higgs potential, one has the following 
relations 



V = m\m\ — mg < 0, 



(4.26) 



C = m\ + m^ — 2|m3| > 0. 



(4.27) 




(4.28) 



(4.29) 
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where tan/3 = ^2/^1, /i^ = rn"^ + = 1,2) and Sj is the loop corrections arising 
from the loop Higgs potential AVi. Taking /x^ and as the input parameters, one 
can solve for tan (3 and Mz 

sin(2/3) = (4.30) 



1 2 _ /xi - /i2 tan^ /? 
2^^^^" tan2/?-l 



M| = ^\7^r; / ■ (4.31) 



One can also treat the tan (3 and Mz as the input parameters, and eliminate two of 
the GUT scale parameters, say Bq and /iq, using Eqs. (14.301 I4.3ip . Since only /i^ 
enters the Eqs. (14. 30[ 14^311) . the sign of fi is undetermined. Thus one has that the low 
energy physics depends on the following parameters 

^0, ^1/2, ^0, tan/5, sign{fi) (4.32) 

where ttlq is the universal scalar mass, mi/2 is the universal gaugino mass, Aq is 
the universal trilinear coupling, tan/5 is the ratio of the two Higgs VEVs in the 
MSSM, and fi is the Higgs mixing parameter that enters via the term fiHiH2 in the 
superpotential. 

4.2 Hyperbolic Branch (HB) of Radiative Symme- 
try Breaking 

The symmetry breaking condition Eq. (14.311) can be rewritten as [l5l US] 

* = i + (4.33) 
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\ 2 2 

$-1 = 4—^ (4.34) 



where the new parameter 

is introduced for the purpose of the measure of naturalness. Using the radiative 
electroweak symmetry breaking constraint and ignoring the 6-quark couphngs, we 
may express the parameter $ as 



1 /moV /AoV /mi/^y /rr^X A/iL 

I^WJ ^'^[wj ^^+Im^J ^^+l^i^J^^+^ 



where 



tan^ i3 -1 
1 



= rrS^^' (4-37) 



^ tan2/5 



— - etan^/?) , (4.38) 



. o Si — Eo tan^ /3 



Here Dq = 1 — {mt/mjY with ~ 200 sin/? GeV, and e, /, (7, k are defined in [i7] . 
To investigate the hmits of mo and mi/2 consistent with the symmetry breaking for 
some given value of $, we rewrite the Eq. (14. 2 p as 

C,ml + C3m;% + C'.Al + A/^L, = M| (^$0 + , (4.41) 



where 



^'1/2 = ^1/2 + ^, C'2 = C2-^^. (4.42) 
For small to moderate values of tan/3, when the loop corrections are small and 
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Ci, , and C3 are all positive from the renormalization group analysis, Eq. (14. 2 p can 
be rewritten as 

^'1/2 ml Al 



2 + (4-43) 

Here one finds that the radiative symmetry breaking demands that the allowed set 
of soft parameters lie on the surface of an Ellipsoid for a fixed value of /i. However, 
for the case with large values of tan/5, the loop corrections to fi become significant. 
In this case, the size of the loop corrections depends sharply on the scale Qq where 
the minimization of the effective potential is carried out. If we choose the scale at 
which the loop corrections are minimized, or even vanish, the loop corrections can 
be omitted for the analysis again except that the sign of Ci can be flipped for some 
region of the parameter space. Typically the scale is not distant from the average of 
the smallest and largest sparticle masses. Under these conditions, the minimization 
condition takes the form 

"^/^ -±1 (4.44) 



where 



. |(^o + l/4)M|-C-Ag^^ ^^^^^^ 



C3 

^2 ^ |(£o + l/4^1^C|A||_ ^^^^g^ 



The set of parameters which satisfy the above relation lie on the surface of a Hyper- 
boloid, and thus this branch is known as the Hyperbolic Branch (HB). The interesting 
thing about the HB region is that rriQ and mi/2 can become quite large with mo lying 
in the multi TeV region even with small fine tuning [15] . 

We note in passing that while the phenomena with soft breaking discussed above 
are within the framework of supergravity models many aspects of these results trans- 
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late to soft breaking within the framework of heterotic string models and for models 
based on intersecting D branes (see, e.g., [l8l H9| [50 | [STl [52] ). 

4.3 Sparticle Masses 



Now let us list the mass matrices for the sparticles in the MSSM. Due to the effect of 
electroweak symmetry breaking, the higgsinos and the electroweak gauginos are mixed 
to form the mass eigenstates called neutralinos and charginos. For the neutralinos, 
the mass matrix takes the form 



Ml 






M2 



V 



-Mzswcp Mzcwcp 
Mzswsp -Mzcwsp 



-Mzswcp Mzswsp 
MzCwCp -MzCwSp 











(4.47) 



where 9^ is the weak angle, sw = sin 6*1^, cw = cos 9w, = sin/3, and = cos (3. 
For the charginos, the mass matrix takes the form 



v^M, 



WC/3 



(4.48) 



The up squark (mass)^ matrix at the electroweak scale is given by 



Mi 



Ml + m„2 + M|(i - Q„s^) cos 2/3 
m„(A„ - /i* cot/?) 



m„(v4; - /icot (3) 
ml + + M|Q„s2^ cos 2(3 



(4.49) 



where Qu = I- And the down squark (mass)^ matrix at the electroweak scale is given 
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by 



a 



+ mrf2 - M|(i + Qaslr) cos 2/3 



V 



md{Ad - 11* tan P) 



md{A* - //tan/3) 
+ ma^ + Mfg^s^^ cos 2/3 



(4.50) 



where Qd — —\- In deducing the above, the relations /iu,d,e = ^,d,e^u,d,e are used for 
the trihnear terms. Finally, the slepton (mass)^ matrix is given by 



Ml 



\ 



m2-M|(|-s2,) cos 2/3 
me{Ae - //*tan/3) 



me{Al — /itan/i) 
mj + ml- Mls"^ cos 2/3 



(4.51) 



4.4 Dark Matter in Supergravity 



In most of the allowed parameter space of the mSUGRA model consistent with ra- 
diative breaking of the EW symmetry, the hghtest neutrahno is the LSP and hence 
a candidate for cold dark matter. We briefly discuss the computation of the relic 
density of neutralino = Px/ Pc where p is the mass density of relic neutralinos in 
the universe and Pc is the critical mass density needed to close the universe, i.e. 



ml 



SttG 



N 



(4.52) 



Here Hq is the Hubble parameter at current time and Gjv is the newtonian constant. 
Numerically, one has 



Pc = 1.9/io X 10 ^^gm/cm 



(4.53) 



where ho is now the Hubble parameter in unit of 100 km/secMpc. The current value 
of ho is ho — 0.7 ± 0.013. In the analysis of il^^o we need to solve the Boltzman 
equation for n, the number density of neutralinos in the early universe, which is given 
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by [53] 



^ = -?,Hn - {av){n^ - nl). (4.54) 



In the above, no is the value of n at thermal equilibrium, {av) is the thermal average 
of the neutralino annihilation cross section o"(x'j'Xi ~^ ^) ^"^^ is the relative Xi 
velocity, and H is the Hubble parameter at time t. In the computation of the thermal 
average one can assume that the neutralinos are non-relativistic and thus one can 
approximate {av) by the relation 

, , {'^ dvv'^{av)e~'"'^ /'^'^ 

Here x is defined by x = kT/m^ where T is the temperature and k is Boltzman 
constant. A solution to the Boltzman equation gives 



n,/4-2.5xlO-"(|)^(4)^U (4.56) 



where is the current microwave background temperature, is the "freeze out" 
temperature corresponding to the temperature where the annihilation rate becomes 
smaller than the expansion rate, so that x\ decouples from the background, is 
typically small with a value x/ ~ 0.04. Nj is the number of degrees of freedom at 
freeze out and typically Nj ~ 289.5/8 [54j. The factor (T^/T^)^ is estimated to be 
(T^/T-x)^ ~ 18.5 (for more recent evaluation, see [55]). Finally, J{xf) is the integral 
defined by 

J{xf)= [ ^ dx{av)x (GeV-2). (4.57) 
Jo 

Away from the poles one can carry out a power series expansion for av so that 
av = a + bv'^/6 + ■■. In this case the thermal average is straightforward. However, 
the above approximation is invalid near poles [56], EH EHl EH] and the integration 
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becomes tricky since one has a double integration over a pole. However, it is possible 
to overcome this problem by an interchange in the order of integration as discussed 
in [581 [59]. 

A realistic computation of the relic density in supersymmetric models is however, 
more complicated as co-annihilations contribute to the relic density [56[ [601 EB [62] 
For instance, one can have co-annihilations involving staus f so that 

fx tZ, rh, r7 (4.58) 
rf* ^ fJi,W+W~,ZZ,jZ,jj (4.59) 
ff TT (4.60) 

fi,{t^T) ^ Tii (4.61) 



For the case of co-annihilation, one must consider the total density n = ^ • rii 
where i runs over all the sparticles that enter in the co-annihilations, where n now 
obeys the equation 

— = -3Hn - {a,sv,,i) {n^ - n^J (4.62) 

where 

(Tcff = ^0"ij7i7j. (4.63) 

Here 0"^ is the cross section of annihilation of particles i and j, and 7j = nl^/rieq 
where n'q refers to the number density of sparticle i at thermal equilibrium. The 
relic density for these processes requires numerical integration programs and we use 
micrOMEGAs [63] in our analysis. We note in passing that relic density is affected 
by Yukawa unification (see, e.g. [64j and the references therein) but we do not take 
such effects into account here. 
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4.5 CP violation 



MSSM contains many sources of CP violation which arise from the soft breaking 
sector of the theory (for a review of CP violation in SUGRA, strings and branes see 
[55]). The number of phases in SUGRA models is reduced drastically. Specifically 
in mSUGRA one has just two CP phases which can be chosen to be the phase of 
the trilinear coupling A and the phase of fi. With non-universalities one can bring 
in new phases. A similar situation holds in string and D brane models. A stringent 
constraint on CP phases arises from the electric dipole moment (EDM) of the electron 
and of the neutron (see, e.g., [55]) which would naively imply that SUSY phases are 
all very small. However, this need not be the case [67] because of the cancellation 
mechanism [68]. With the cancellation mechanism, one finds that the phases can be 
large and at the same time one can satisfy the EDM constraints. Such phases can 
affect low energy phenomena such as the the gaugino and sfermion masses and Higgs 
masses [69] . CP phases can affect LHC signatures. However, in this study we do not 
take the effect of CP phases into account. 



4.6 Proton stability 

Another constraint on unified models of particle interactions arises from proton decay 
constraints (for a recent review see [TOj). In grand unified models and also in string 
and brane models, one has several sources of proton decay. First of all, both in super- 
symmetric as well as in non-supersymmetric grand unified models one has baryon and 
lepton number violating dimension six operators due to the exchange of lepto-quarks. 
The most prominent decay mode from these is p — > e^n^ and the current limit on this 
decay is r(p —>■ e+vr'^) > 1.6 x 10'^^ yr. In supersymmetric theories the most dominant 
decay mode is p — »• i/K'^ and the current limit on it is r(p vK^) > 2.3 x 10^^ yr. 
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This latter limit, i.e., on the mode vK^ ^ puts very stringent limits on grand unified 
models and on string models. At the same time these modes are subject to a much 
greater degree of model dependence because of the unknown nature of physics at high 
scales. In contrast the soft parameters are known to be independent of the high scale, 
specifically of the GUT scale |6]. For this reason we will not consider any specific 
high scale model but rather focus on weak scale supersymmetry which is determined 
by the soft parameters independent of the GUT scale [6]. 
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Chapter 5 

Experimental Constraints on 
Unified Models 

Below we summarize the relevant constraints from collider and from astrophysical 
data that are applied throughout the analysis unless stated otherwise. 

1. WMAP 3 year data: The lightest R-Parity odd supersymmetric particle (LSP) 
is assumed charge neutral. The constraint on the relic abundance of dark matter 
under the assumption that the relic abundance of neutralinos is the dominant 
component places the bound: 0.0855 < ^^oJi^ < 0.1189 {2a) [71] . 

2. As is well known sparticle loop exchanges make a contribution to the FCNC 
process b s'y which is of the same order as the Standard Model contributions 
(for an update of SUSY contributions see [72]). The experimental limits on 
b s'-f impose severe constraints on the SUSY parameter space and we use 
here the constraints from the Heavy Flavor Averaging Group (HFAG) [73] along 
with the BABAR, Belle and CLEO experimental results: Br{b 57) = (355 ± 
24l?o ± 3) X 10-^ A new estimate of Br{B X^-f) at 0{a^) gives [71] Br{b 
57) = (3.15±0.23) X 10"^ which moves the previous SM mean value of 3.6 x 10~^ 
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a bit lower. In order to accommodate this recent analysis on the SM mean, as 
well as the previous analysis, we have taken a wider 3.5a error corridor around 
the HFAG value in our numerical analysis. The total Br{B — > Xg"^) including 
the sum of SM and SUSY contributions are constrained by this corridor. With a 
2cr corridor, while some of the allowed points in our analysis will be eliminated, 
the main results of our pattern analysis remain unchanged. 

3. The process Bs — > can become significant for large tan/5 since the decay 
has a leading tan^ P [75] dependence and thus large tan (3 could be constrained 
by the experimental limit Br{Bs < 1.5 x 10"^ (90% CL), 2.0 x 
10^^ (95% CL) |136j . This limit has just recently been updated [TB] and gives 
Br{Bs — s> < 1.2 X 10^'' (95% CL). Preliminary analyses [77] have reported 
the possibility of even more stringent constraints by a factor of 10. We take a 
more conservative approach in this analysis and allow model points subject to 
the bound Br{Bs IJ'^I^ ) < 9 x 10^^ (for a review see [78]). 

4. Additionally, we also impose a lower limit on the lightest CP even Higgs boson 
mass. For the Standard Model like Higgs boson this limit is ~ 114.4 GeV [79], 
while a limit of 108.2 GeV at 95% CL is set on the production of an invisibly 
decaying Standard Model like Higgs by OPAL For the MSSM we take the 
constraint to be nih > 100 GeV. A relaxation of the light Higgs mass constraint 
by 8 - 10 GeV affects mainly the analysis of SUGRA models where the stop 
mass can be light. However, light stops are possible even with the strictest 
imposition of the LEP bounds on the SM Higgs Boson. We take the other 
sparticle mass constraints to be m~± > 104.5 GeV [81] for the lighter chargino, 
m^^ > 101.5 GeV for the lighter stop, and m^r^ > 98.8 GeV for the lighter stau. 
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In addition to the above one may also consider the constraints from the anomalous 
magnetic moment of the muon. It is known that the supersymmetric electroweak 
corrections to — 2 can be as large or larger than the Standard Model electroweak 
corrections [82]. The implications of recent experimental data has been discussed 
in several works (see, e.g. [83]). As in [33], here we use a rather conservative bound 
-11.4 X 10-10 < - 2 < 9.4 X 10-9. 
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Chapter 6 

The SUGRA Sparticle Patterns 



In this chapter, we discuss the possible sparticle mass hierarchical patterns arising 
from SUGRA models utilizing the Monte Carlo simulations. We also discuss the 
correlation between the sparticle patterns and the nature of the soft parameter space. 
Finally we give a collection of benchmarks for each of the sparticle patterns analyzed 
here. 



6.1 The Sparticle Landscape 

The analysis proceeds by specifying the model input parameters at the GUT scale, 
Mg ~ 2 X 10^^ GeV, (for our analysis, no flavor mixing is considered at the GUT scale) 
and using the renormalization group equations (RGBs) to predict the sparticle masses 
and mixing angles at the electroweak scale. The RGE code used to obtain the mass 
spectrum is SuSpect 2.34 [H], which is the default RGE calculator in MicrOMEGAs 
version 2.0.7 [63]. We have also investigated other RGE programs including ISAS- 
UGRA/ISAJET [85], SPheno [86] and SOFTSUSY [87]. We have cross checked our 
analysis using different codes and find no significant disagreement in most regions of 
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the parameter space. The largest sensitivity appears to arise for the case of large tan (3 
and the analysis is also quite sensitive to the running bottom mass and to the top 
pole mass (we take m^^(m;,) = 4.23 GeV and m((pole) = 170.9 GeV in this analysis). 
Such sensitivities and their implications for the analysis of relic density calculations 
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Table 6.1: Hierarchical mass patterns for the four lightest sparticles in mSUGRA 
when /i < and /i > 0. The patterns can be classified according to the next to the 
lightest sparticle. For the mSUGRA analysis the next to the lightest sparticle is found 
to be either a chargino, a stau, a stop, a CP even/odd Higgs, or the next lightest 
neutralino X2- The notation A/H stands for either A or H. In mSP14-mSP16 it is 
possible that the Higgses become lighter than the LSP. Y stands for appearance of 
the pattern for the sub case. 



are well known in the literature |SH] and a detailed comparison for various codes can 
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be found in Refs. ([89], [90], [91], [92]). 



6.1.1 The mSUGRA landscape for the 4 hghtest sparticles 

As discussed in chapter one mSUGRA model is a point in a 4 dimensional pa- 
rameter space spanned by the soft parameters mo, mi/2, Aq, tan/5, and the sign of 
fi. Typically scans of the parameter space are done by taking a vanishing trilinear 
coupling, and/or by looking at fixed values of tan/5 while varying (mo, mi/2). In this 
work we carry out a random scan in the 4-D input parameter space for both signs 
of fi with Monte Carlo simulations using flat priors under the following ranges of the 
input parameters 



0<mo<4TeV, < mi/2 < 2 TeV |Ao/mo| < 10, 1 < tan/3 < 60. (6.1) 



Since SUGRA models with /i > are favored by the experimental constraints much of 
the analysis presented here focuses on this case. Specifically for the /i > mSUGRA 
case, we perform a scan of the parameter space with a total of 2 x 10^ trial parameter 
points. We delineate the patterns that emerge for the first four lightest sparticles. 
Here we find that at least sixteen hierarchical mass patterns emerge which are la- 
beled as mSPs (minimal SUGRA Pattern). These mSPs can be generally classified 
according to the type of particle which is next heavier than the LSP, and we find 
four classes of patterns in mSUGRA: the chargino patterns (CP), the stau patterns 
(SUP), the stop patterns (SOP), and the Higgs patterns (HP), as exhibited below 

1. Chargino patterns (CP) : mSPl, mSP2, mSP3, mSP4; 

2. Stau patterns (SUP) : mSP5, mSP6, mSP7, mSP8, mSP9, mSPlO; 

3. Stop patterns (SOP) : mSPll, mSP12, mSP13; 
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mSUGRA Pattern (mSP) 

Figure 6.1: Distribution of the surviving hierarchical mass patterns in the landscape 
for the mSUGRA model with /i > (light) and fi < (dark), under various constraints 
as discussed in the text. 

4. Higgs patterns (HP) : mSP14, mSP15, mSP16. 

The hierarchical mass patterns mSPl-mSP16 are defined in Table (16. ip . We note 
that the pattern mSP7 appears in the analyses of [931 EH ES]. We also performed 
a similar scan for the mSUGRA with /x < case using the Monte Carlo simulation 
with flat priors and the same parameter ranges as specified in Eq. (16.11) . Most of the 
mSP patterns that appear in the fi > case also appear in the fi < case (see Table 
(16. ip ). However, in addition one finds new patterns shown below 

1. Stau patterns (SUP) : mSP17, mSP18, mSP19; 

2. Stop patterns (SOP) : mSP20, mSP21; 

3. Neutralino patterns (NP) : mSP22. 

We note that the analysis of Ref. [96] has a sparticle spectrum which corresponds 
to mSPll and contains light stops. Light stops have also been discussed recently in 

45 



Snowmass 



SPSla, SPSlb, SPSS 



SPS2 



SPSS 



SPS4, SPS6 



mSP 



mSP7 



mSPl 



mSP5 



mSP3 



Post-WMAP3 



A', B', C, D', G', H', J', M' 



I',L' 



E' 



K' 



mSP 



mSP5 



mSP7 



mSPl 



mSP6 





CMS LM/HM 


mSP 






LMl, LM6, HMl 


mSP5 






LM2, LM5, HM2 


mSP7 






LM3, LM7, LM8, LM9, LMIO, HM4 


mSPl 






LM4, HM3 


mSP3 





Table 6.2: Mapping between the mSPs and the Snowmass, Post-WMAP3, and CMS 
benchmark points. The points B' = LMl, /' = LM2, C = LM6. HMl in SuSpect has 
rriy^o > m^f^, but this is not the case for ISAJET, SPheno, and SOFTSUSY. Among the 
CMS benchmarks, only LMl, LM2, LM6, and HMl, HM2 are capable of giving the 
correct relic density. Thus the mapping above applies only to the mass pattern, while 
all of our mSP and NUSP benchmark points satisfy the relic density constraints from 
MicrOMEGAs with SuSpect. The CMS test points do a better job of representing 
mSPl which is the dominant pattern found in our analysis. There are no HP test 
points or SOP test points in any of the previous works. 



While the earlier works which advocated benchmark points and slopes made good 
progress in systematizing the search for supersymmetry, we find that they do not cover 
the more broad set of possible mass hierarchies we discuss here. That is, many of the 
mSP patterns do not appear in the earlier works that advocated benchmark points 
for SUSY searches. For example, the Snowmass mSUGRA points (labeled SPS) [98] 
and the Post-WMAP benchmark points of [99], make up only a small fraction of the 
possible mass hierarchies listed in Table (16. ip . The CMS benchmarks classified as 
Low Mass (LM) and High Mass (HM) [IDQ] (for a recent review see [Ml EDS]) does 
a good job covering the mSPl pattern which appears as the most dominant pattern 
in our analysis, but there are no Higgs patterns or stop patterns discussed in the 
CMS benchmarks as well as in SPS or in Post-WMAP benchmarks. We exhibit the 
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mapping of mSPs with other benchmarks points in a tabular form in Table (16.20 . 

In Fig. (16. ip we give the relative distribution of these hierarchies found in our 
Monte Carlo scan. Because the scan is done randomly within the soft parameter 
space, the distribution of sparticle patterns in Fig. (16.11) represents the probability of 
finding these patterns in the parameter space. The most common patterns found here 
are CPs and SUPs, especially mSPl and mSP5. However there exists a significant 
region of the parameter space where SOPs and HPs can be realized. The percentages 
of occurrence of the various patterns in the mSUGRA landscape for both /i positive 
and fi negative are exhibited in Fig. (16.11) . The analysis of Fig. (16.11) shows that 
the chargino patterns (CP) are the most dominant patterns, followed by the stau 
patterns (SUP), the stop patterns (SOP), and the Higgs patterns (HP). In contrast, 
most emphasis in the literature, specifically in the context of relic density analysis, has 
focused on the stau patterns, with much less attention on other patterns. Specifically 
the Higgs patterns have hardly been investigated or discussed. The exceptions to this, 
in the context of the Higgs patterns, are the more recent works of Refs. [221 123] , and 
similar mass ranges for the Higgs bosons have been studied in |103j (see also |104] ). 

6.1.2 The landscape of the 4 Hghtest sparticles in NUSUGRA 

Next we discuss the landscape of the 4 lightest sparticles for the case of nonuni- 
versal supergravity models. Here we consider nonuniversalities in the Higgs sector 
(NUH), in the third generation sector (NU3), and in the gaugino sector (NUG). Such 
nonuniversalities appear quite naturally in supergravity models with a non-minimal 
Kahler potential, and in string and D-Brane models. The parametrization of the 
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nonuniversalities is given by 



NUH : Mh„ =mo(l + 5i^J, Mh, = m^il + 5h,), 

NU3 : M,3 = mo(l + 5g3), M^za^ = m^{l + 5tm), (6-2) 

NUG : Mi = mi/2, Ma.s = mi/2(l + ^a/^.s)- 

In the above 5h^ and 5h^ define the nonuniversahties for the up and down Higgs 
mass parameters, M^s is the left-handed squark mass for the 3rd generation, and M^s 
(Mrfs) are the right-handed u-squark (d-squark) masses for the 3rd generation. The 
nonuniversahties in the gaugino sector are parameterized here by and 5ms- We 
have carried out a Monte Carlo scan with flat priors using 10^ model points in each 
of the three types of NUSUGRA models, taking the same input parameter ranges as 
specified in Eq. (16. ip and —0.9 ^5^1. Almost all of the mSP patterns seen for the 
mSUGRA cases were found in supergravity models with nonuniversal soft breaking, 
as the mSUGRA model is contained within the nonuniversal supergravity models. In 
addition we find many new patterns labeled NUSPs (nonuniversal SUGRA pattern), 
and they are exhibited in Table (16.31) . As in the mSUGRA case one finds several 
pattern classes, CPs, SUPs, SOPs, and HPs as exhibited below. In addition, we find 
several Gluino patterns (CP) where the gluino is the NLSP. 

1. Chargino patterns (CP) : NUSPl, NUSP2, NUSP3, NUSP4; 

2. Stau patterns (SUP) : NUSPS, NUSP6, NUSP7, NUSPS, NUSP9; 

3. Stop patterns (SOP) : NUSPIO, NUSPll; 

4. Higgs patterns (HP) : NUSP12; 

5. Gluino patterns (CP) : NUSP13, NUSP14, NUSP15. 



48 



It is interesting to note that for the 4 sparticle landscape we find saturation in the 
number of mass hierarchies that are present. For example, for the case > in 
mSUGRA , increasing the soft parameter scan from 1 x 10^ parameter model points 
to 2 X 10^ model points does not increase the number of 4 sparticle patterns. In 
this context it becomes relevant to examine as to what degree the relic density and 
other experimental constraints play a role in constraining the parameter space and 
thus reducing the number of patterns. This is exhibited in Table (16. 4p where we 
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Table 6.3: New 4 sparticle mass patterns that arise in NUSUGRA over and above 
the mSP patterns of Table (16. ip . These are labeled nonuniversal SUGRA patterns 
(NUSP) and at least 15 new patterns are seen to emerge which are denoted by NUSPl- 
NUSP15. 



demonstrate how the relic density and the other experimental constraints decrease the 
number of admissible model points in the allowed parameter space for the mSUGRA 
models with both /i > and /i < 0, and also for the cases with nonuniversalities in 
the Higgs sector, nonuniversalities in the third generation sector, and with nonuniver- 
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salities in the gaugino sector. In each case we start with 10^ model points at the GUT 
scale, and find that the electroweak symmetry breaking constraints reduce the num- 
ber of viable models to about 1/4 of what we started with. We find that the allowed 
number of models translates into SUGRA mass patterns which are typically less than 
100. The admissible set of parameter points reduces drastically when the rehc density 
constraints are imposed and are then found to typically reduce the number of models 
by a factor of about 200 or more, with a reduction in the number of allowed pat- 
terns by a factor of 2 or more. Inclusion of all other experimental constraints further 
reduces the number of admissible points by a factor between 30% and 50%, with a 
corresponding reduction in the number of patterns by up to 40%. The above anal- 
ysis shows that there is an enormous reduction in the number of admissible models 
and the corresponding number of hierarchical mass patterns after the constraints of 
radiative breaking of the electroweak symmetry, relic density constraints, and other 
experimental constraints are imposed. 



Model 


Trial 


Output 


No. of 


Relic Density 


No. of 


All 


No. of 


Type 


Models 


Models 


Patterns 


Constraints 


Patterns 


Constraints 


Patterns 


mSUGRA(^ > 0) 


10'' 


265,875 


55 


1,360 


22 


902 


16 


mSUGRA(At < 0) 


10'' 


226,991 


63 


1,000 


31 


487 


18 


NUH(/i > 0) 


10" 


222,023 


59 


1,024 


24 


724 


15 


NU3(/t > 0) 


10" 


229,928 


73 


970 


28 


650 


20 


NUG(/i > 0) 


10" 


273,846 


103 


1,788 


36 


1,294 


28 



Table 6.4: An analysis of mass patterns for the four lightest sparticles. Exhibited 
here are the model type, the number of trial input points for each model, the number 
surviving the radiative electroweak symmetry breaking scheme as given by SuSpect 
(column 3) , the number surviving when the relic density constraints are applied with 
MicrOMEGAs (column 5), the number surviving with inclusion of all experimental 
collider constraints (column 7) , along with the corresponding number of hierarchical 
mass patterns in each case (column 8). 
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6.1.3 Hierarchical patterns for the full sparticle spectrum 

We discuss now the number of hierarchical mass patterns for the full set of 32 sparti- 
cles in SUGRA models when the constraints of electroweak symmetry, relic density, 
and other experimental constraints are imposed. The result of the analysis is given 
in Fig. (16. 2p and Table (16. 5p . The left panel of Fig. (16. 2p shows the number of hi- 




Number of mSUGRA Models Number of mSUGRA Models 

Figure 6.2: Saturation of Sparticle patterns. 



Models [No.] 


No. after constraints 


No. of patterns 


mSUGRA(/^ > 0) [10''] 


902 


505 


mSUGRA(/^ < 0) [10''] 


487 


268 


NUH(/i > 0) [10'=] 


724 


517 


NU3(Ai > 0) [10'^] 


650 


528 


NUG(At > 0) [10"] 


1294 


1092 


All Above[5 x 10"] 


4057 


2557 



Table 6.5: Exhibition of mass patterns and models with various constraints. Column 
1 shows one million input parameter points for each of the models investigated, and 
the number surviving all the constraints are exhibited in column 2, while column 3 
gives the number of hierarchical patterns. 



erarchical mass patterns for 32 sparticles vs the number of trial points for mSUGRA 
models which survive the electroweak symmetry breaking constraints, the relic den- 
sity and all other experimental constraints. The number of hierarchical mass patterns 
show a trend towards saturation. In the right panel, a similar phenomenon is seen in 
the ratio between the number of patterns over the number of surviving trial points 
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in mSUGRA models. Here one finds that increasing the number of model points in 
the scan does increase the number of patterns. However, the ratio of the number of 
patterns to the total number of models that survive all the constraints from the scan 
decreases sharply as shown in the right panel of Fig. (16. 2p . This means that although 
saturation is not yet achieved one is moving fast towards achieving saturation with a 
relatively small number of allowed patterns for all the 32 sparticles within SUGRA 
models consistent with the various experimental constraints. The analysis of Table 
(16. 5p shows that the number of allowed patterns for the 32 sparticles, which in the 
MSSM without the SUGRA framework can be as large as 0(10^^) or larger, reduces 
rather drastically when various constraints are applied in supergravity models. 

The Table (16.51) exhibits a dramatic reduction of the landscape from upward of 
~ 0(10^^) hierarchical mass patterns for the 32 sparticle masses to a much smaller 
number when the electroweak symmetry breaking constraints, the relic density con- 
straints, and other experimental constraints are applied. We note that some patterns 
are repeated as we move across different model types listed in the first column of Ta- 
ble (16.51) . Thus the total number of patterns listed at the bottom of the last column 
of this table is smaller than the sum of patterns listed above in that column. We 
note that the precise number and nature of the patterns are dependent on the input 
parameters such as the top mass and a significant shift in the input values could 
modify the pattern structure. 

6.2 Sparticle Patterns & Nature of Soft Breaking 

It is interesting to ask if the patterns can be traced back to some specific regions of 
the parameter of soft breaking from where they originate. This indeed is the case, 
at least, for some of the patterns. The analysis illustrating the origin of the patterns 
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in the parameter space is given in Fig. (16.31) . The dispersion of mSPs arising in 
mSUGRA in the tan (3 vs Ao/mo plane (left panels), and in the mo vs mi/2 plane 
(right panels) for the /i > case (upper panels) and /i < case (lower panels). The 
analysis is based on a scan of 10^ trial model points with flat priors in the ranges 
mo < 4 TeV, mi/2 < 2 TeV, 1 < tan/3 < 60, and |y4o/mo| < 10. mSPl is confined 
to the region where |Ao/mo| < 2. For the case /i < 0, no HPs are seen, and also, no 
model points survive in the region where tan/5 > 50 in contrast to the ^ > case 
where there is a significant number for tan f3 > 45. Many interesting observations can 
be made from these spectral decompositions. For example, a significant set of the 
mSPl (CP) models lie in the region l^o/mol < 2 and correspond to the Hyperbolic 
branch/Focus Point (HB/FP) |15] regions, while most of the SOPs have a rather large 
ratio of Ao/mQ with the satisfaction of REWSB. In this analysis we require that there 
be no charge or color breaking (CCB) |105[ I106j at the electroweak scale. We note in 
passing that it has been argued that even if the true minimum is not color or charge 
preserving, the early universe is likely to occupy the CCB preserving minimum and 
such minima may still be acceptable if the tunneling lifetime from the false to the 
true vacuum is much greater than the present age of the universe |107j . Next, we note 
that for the mSUGRA /i > case, the region around tan j3 = 50 has a large number 
of models that can be realized, while the region around tan f3 = 30 has far less model 
points. We also note that most of the HPs reside only in the very high tan j3 region 
in mSUGRA, but this situation can be changed significantly in the NUH case where 
HP points can be realized in the tan/3 region as low as tan/3 ~ 20. In the mo vs mi/2 
plane, one finds that most of CPs and HPs have a larger universal scalar mass than 
most of the SUPs and SOPs. 

Often in the literature one limits the analysis by fixing specific values of Aq and 
tan/3. For Aq the value most investigated is Aq = 0. However, constraining the 
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Figure 6.3: Spectrum decomposition of soft parameter space. Exhibited are the 
landscape of sparticle mass spectra in the planes (I) tan/3 vs Ao/niQ and (II) mo vs 
mi/2, when the soft parameters are allowed to vary in the ranges given in Eq. (16.11) . 



values of Aq or of tan/3 artificially eliminates a very significant part of the allowed 
parameter space where all the relevant constraints (the REWSB constraint as well 
as the relic density and the experimental constraints) can be satisfied as seen in 
Fig. (16.31) . One can extract the familiar plots one finds in the literature where Aq 
and tan/? are constrained from a reduction of the top-right panel of Fig. (16. 3p . The 
results of this reduction are shown in Fig. (16.41) with a focused scan in specific regions 
of the soft parameter space. Fig. (16.41) shows a dispersion of patterns in the mo vs 
mi/2 plane for fixed values of tan/3 and Ao/niQ. The region scanned is in the range 
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Figure 6.4: Dispersion of patterns in some special 2D surface of mSUGRA space. 

mo < 4 TeV and mi/2 < 2 TeV with a 10 GeV increment for each mass. Only a subset 
of the allowed parameter points relative to Fig. (16. 3p remain, since the scans are on 
constrained surfaces in the mSUGRA parameter space. Specifically the bottom-left 
and top-right panels of Fig. (16.41) show the familiar stau coannihilation [lOSj I109[ 193] 
regions and the HB/FP branch, the bottom- right panel gives the stau coannihilation 
region and the stop coannihilation region because of the relatively large Aq value, and 
the top-left panel is of the form seen in the works of Djouadi et al. [33j where the 
Higgs funnel plays an important role in the satisfaction of the relic density. 

A similar analysis for the nonuniversal case is given in Fig. (16.51) which exhibits 
the NUSPs and mSPs for the NUH, NU3, and NUG models in the tan/5 vs Ao/ruo 
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Figure 6.5: Exhibition of NUSPs and mSPs arising from NUSUGRA models. 



plane. The range of SUGRA parameters are the same as the case mSUGRA (/i > 0). 
One may notice that the mSPl points arising from NU models lie in a relatively 
larger Ao/mQ region. Most of the models in NU cases are still mSPs, and among the 
NUSPs, only two patterns have a relatively large population, these being NUSPl and 
NUSP13. One may also notice that in NUH case, the HPs can exist in a low tan/5 
region as opposed to the mSUGRA case where HPs can either exist in the large tan f3 
region (yU > 0) or are totally eliminated < 0). Among the NUSPs the dominant 
patterns are NUSPl (CP) and NUSP13 (GP), which are seen to arise from the model 
with nonuniversalities in the gaugino sector, i.e., the NUG model. 
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6.3 Benchmarks for Sparticle Patterns 

As discussed already, many of the sparticle mass patterns discussed in this analysis do 
not appear in the Snowmass, Post-WMAP, and CMS benchmark points. With some 
of these mSP and NUSP having a significant probability of occurrence, we therefore 
provide a larger set of benchmark points for the various patterns in different SUGRA 
scenarios. In Table (16.61) . we give one benchmark point for each mSP pattern that 
are discovered in mSUGRA with /i > 0. 



mSUGRA 


mo 




Ao 


tan/3 


KQ) 


LSP 


LCP 


Pattern 


(GeV) 


(GeV) 


(GeV) 


Vu/Vd 
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(GeV) 


mSPl: 
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mSP2: 
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mSP8: 




1880 
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4075 


54.8 


1141 


373.1 


379.6 


mSP9: 




667 


1154 


-125 


51 


1257 


499.2 


501.8 


mSPlO 




336 


772 


-3074 


10.8 


1695 


329.2 


331.7 


mSPll 




871 


1031 


-4355 


10 


2306 


447.1 


491.5 


mSP12 
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1671 


-6855 


10 


3593 


741.2 


791.8 


mSP13 




524 


800 


-3315 


15 


1782 


342.7 


383.8 


mSP14 




1036 


562 


500 


53.5 


560 


236.2 


399.1 


mSP15 




1113 


758 


1097 


51.6 


724 


321.1 


595.9 


mSP16 




525 


450 


641 


56 


484 


184.6 


257.9 



Table 6.6: Benchmarks using SUSPECT 2.3 with one point for each mass pat- 
tern mSPl-mSP16. Also given are the neutralino LSP (Lightest SUSY (R parity 
odd) Particle), and the Lightest Charged Particle (LCP) masses. We take /i > 0, 
rrib^imb) = 4.23 GeV, a,^(Mz) = .1172, and mt(pole) = 170.9 GeV. At least five 
LCP from these benchmarks will be accessible at the International Linear Collider 
(ILC). 



We also provide a large collection of benchmark points which are exhibited in 
Tables (115.11 115.21 115.31 115.41 115. 5p in the Appendix. Each of these benchmarks 
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satisfies the relic density and other experimental constraints with SuSpect linked to 
MicrOMEGAs. We have explicitly checked that the first mSP benchmark point in 
each of the tables can be reproduced by using SPheno, and SOFTSUSY by allowing 
minor variations on the input parameters. The benchmarks are chosen to cover wide 
parts of the SUGRA parameter space. We give these benchmarks, several for each 
mass pattern, as the search for SUSY from the point of view of mass patterns has 
important consequences for LHC experimental searches. Some of the patterns are 
correlated with certain well investigated phenomena such as the HB/FP branches 
of REWSB and the stau-neutralino co-annihilation regions. However, many of the 
patterns arise from multiple annihilation processes. 
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Chapter 7 

Sparticle Signatures at the Large 
Hadron Collider 

In this chapter, first we will give the detailed description of our LHC simulations. 
And then we will give our post trigger cuts and all the signatures investigated. 

7.1 LHC Simulations 

After the imposition of all the constraints mentioned in the previous sections, such 
as the relic density constraints from WMAP data, the constraints on the FCNCs, as 
well as mass limits on the sparticle spectrum, we are left with the candidate model 
points for the signature analysis. For each of these model points, a SUSY Les Houches 
Accord (SLHA) file [HO] is interfaced to PYTHIA 6.4.11 [m] through PGS4 [TT2] 
for the computation of SUSY production cross sections and branching fractions. In 
this analysis, for signals, we have generated all of PYTHIA's 2 — 2 SUSY production 
modes using MSEL = 39. More specifically this choice generates 91 SUSY production 
modes including gaugino, squark, slepton, and SUSY Higgs pair production but leaves 
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out singly produced Higgs production. For further details, see A treatment of 

singly produced Higgs production in the context of sparticle mass hierarchies was 
included in the analysis of Ref . [23] . Leading order cross sections from P YTHIA and 
leading order cross sections from PROSPINO 2.0 |113j were cross checked against one 
another for consistency over several regions of the soft parameter space. TAUOLA 
[114] is called by PGS4 for the calculation of tau branching fractions as controlled in 
the PYTHIA parameter card (.pyt) file. 

With PGS4 we use the Level 1 (LI) triggers based on the Compact Muon Solenoid 
detector (CMS) specifications |115[ 1100] and the LHC detector card. Muon isolation 
is controlled by employing the cleaning script in PGS4. We take the experimental 
nomenclature of lepton being defined only as electron or muon and thus distinguish 
electrons and muons from tau leptons. SM backgrounds have been generated with 
QCD multi-jet production due to light quark flavors, heavy flavor jets {bb, tt), Drell- 
Yan, single Z/W production in association with quarks and gluons {Z/W+ jets), and 
ZZ, WZ, WW pair production resulting in multi-leptonic backgrounds. Extraction of 
final state particles from the PGS4 event record is accomplished with a code SMART 
(= SUSY Matrix Routine) written by us [22] which provides an optimized processing 
of PGS4 event data files. The standard criteria for the discovery limit of new signals 
is that the SUSY signals should exceed either Si/iVsM or 10 whichever is larger, i.e., 
NsusY > Max {Sv'iVsM, 10} and such a criteria is imposed where relevant. We have 
also cross checked various results of our analysis with three CMS notes |116l 11171 1118] 
and we have found agreement with these works using SMART and PGS4 for signal 
and backgrounds. 

We note that several works where sparticle signatures are discussed have appeared 
recently |119ill20],ll2H[T22|ll23[[T2l] . However, the issue of hierarchical mass patterns 
and the correlation of signatures with such patterns has not been discussed which is 
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what the analysis of this work investigates. 
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ITlb 


1 r + 1 b-jet 


2Llb 


2 Leptons + 1 b-jet 


2Tlb 


2 r + 1 b-jet 


0L2b 


Lepton + 2 b-jets 


0T2b 


r + 2 b-jets 


lL2b 


1 Lepton + 2 b-jets 


lT2b 


1 r + 2 b-jets 


2L2b 


2 Leptons + 2 b-jets 


2T2b 


2 r + 2 b-jets 


ep 


in IL 


em 


e~ in IL 


mp 


in IL 


mm 


\x~ in IL 


tp 


r+ in IT 


tm 


T~ in IT 


OS 


Opposite Sign Di-Leptons 


Ob 


b-jet 


SS 


Same Sign Di-Leptons 


lb 


1 b-jet 


OSSF 


Opposite Sign Same Flavor Di-Leptons 


2b 


2 b-jets 


SSSF 


Same Sign Same Flavor Di-Leptons 


3b 


3 b-jets 


OST 


Opposite Sign Di-r 


4b 


4 b-jets and more 


SST 


Same Sign Di-r 


TL 


1 r plus 1 Lepton 



Table 7.1: The table gives a list of 40 counting signatures for each point in the SUGRA 
model parameter space. L = e, signifies only electrons and muons. 



7.2 Post Trigger Level Cuts and LHC Signatures 

Generally speaking, there are two kinds of LHC signatures: (i) event counting sig- 
natures, and (ii) kinematical signatures. We have investigated both of these for the 
purpose of discriminating the sparticle mass patterns. We list our event counting 
signatures in Table (17. ip , where we have carried out analyses of a large set of lepton 
-|- jet signals. In our counting procedure, only electron and muon are counted as 
leptons, while tau jets are counted independently. For clarity, from here on, our use 
of 'jet(s)' will exclude tau jets. Thus, for jet identification, we divide jets into two 
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Kinematical signatures 


1. 




2. 


Effective Mass = P^''^^ + T.i 


3. 


Invariant Mass of all jets 


4. 


Invariant Mass of e^e^ pair 


5. 


Invariant Mass of fJ-^fJ-^ pair 


6. 


Invariant Mass of r"*"r~ pair 



Table 7.2: The table give a list kinematical signatures analyzed. 

categories: b-tagged jets and jets without b-tagging, which we simply label as b-jets 
and non-b-jets (see also [119]). There are some counting signatures that only concern 
one class of measurable events, for example, the number of events containing one 
tagged b-jet and any other final state particles. There are also types of signatures of 
final state particles with combinations of two or three different species. For instance, 
one such example would be the number of events in which there is a single lepton and 
a single tau. 

When performing the analysis of event counting, for each SUGRA model point, 
we impose global post trigger cuts to analyze most of our PGS4 data. Below we 
give our default post trigger cuts which are used throughout this thesis unless stated 
otherwise. 

1. In an event, we only select photons, electrons, and muons that have transverse 
momentum > 10 GeV and |?7^| < 2.4, p = (7, e, /i). 

2. Taus which satisfy > 10 GeV and \ri'^\ < 2.0 are selected. 

3. For hadronic jets, only those satisfying P^ > 60 GeV and \r]^ < 3 are selected. 

4. We require a large amount of missing transverse momentum, P™'^^ > 200 GeV. 

5. There are at least two jets that satisfy the Pt and rj cuts. 
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Our default post trigger level cuts are standard and are designed to suppress the 
Standard Model background, and highlight the SUSY events over a broad class of 
models. 

The different kinematical signatures we investigated for the purpose of discrimi- 
nating among sparticle mass patterns are exhibited in Table (17.21) . One may further 
divide the kinematical signatures into two classes: namely those involving transverse 
momentum Pt and those which involve invariant mass. For those involving Py, we 
have investigated missing Pt distributions and the effective mass, the latter being 
the sum of missing Pt and Pt of all jets contained within an event. For the kinemat- 
ical variables using invariant mass, we reconstruct such quantities for four different 
cases, i.e., the invariant mass for all jets, for e^e~ pair, for fi~^fi~ pair, and for r+r~ 
pair. The reconstruction of the invariant mass of r+r^ pair is based on hadronically 
decaying taus (for recent analyses see [93]). 
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Chapter 8 

Event Counting Signatures for 
Sparticles at the LHC 



8.1 Discrimination among mSPs in mSUGRA 

We turn now to a discussion of how one may distinguish among different patterns. 
The analysis begins by considering the 902 model points that survive our mSUGRA 
scan with 10® trial points, and simulating their LHC signals with PGS4 using, for 
illustration, 10 fb~^ of integrated luminosity at the LHC. In our analysis we will 
focus mostly on the counting signatures. Here the most useful counting signature is 
the total number of SUSY events after trigger level cuts and post trigger level cuts 
are imposed. All other counting signatures are normalized with respect to the total 
number of SUSY events passing the cuts and thus appear as fractions lying between 
(0,1) in our figures. To keep the analysis statistically significant, we admit only those 
points in the parameter space that generate at least 500 total SUSY events. 

We give now the details of the analysis. In Fig. (18. II) . we investigate the signature 
space spanned by a variety of signature channels. Top Left: An exhibition of the 
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Figure 8.1: Exhibitions of the mSPs in mSUGRA with /i > in various signature 
channels. 

mSPs in the IL vs OL where the fraction of events to the total number of events in 
each case is plotted. The analysis shows that the Stop Patterns (SOP) appearing on 
the right-bottom corner are easily distinguished from other patterns. The analysis 
shows that SOP has few lepton signals. Top Right and Bottom Left: Plots in the 
signature space with fraction of events with lb vs Ob and 2b vs lb exhibiting the 
separation of CPs and HPs from SOPs and SUPs, with CPs and HPs occupying one 
region, and SOPs and SUPs occupy another in this signature space except for a very 
small overlap. Bottom Right: An exhibition of the mSPs in the signature space with 
the average missing Pt for each parameter point in the mSUGRA parameter space 
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Figure 8.2: An exhibition of how the mSPs can be discriminated within CPs, SUPs. 

along the y-axis and the fraction of events with 06 along the x-axis. The plot shows 
a separation of the CPs and HPs from SOPs and SUPs. Further, mSP4 appears 
isolated in this plot. Most of the CPs and HPs have less than 60% events without 
b-jet content. The ratios for the SUSY models refer to the SUSY signal only. The 
SM point is purely background. The top left panel gives a plot with one signature 
consisting of events with one lepton and the second signature consisting of event with 
no leptons. It is seen that the stop patterns (SOPs) that survive the cuts are confined 
in a small region at the right-bottom corner and have a significant separation from 
all other mSPs. The panel illustrates the negligible leptonic content in stop decays. 
The top-right panel is a plot between two signatures where one signature contains a 
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Figure 8.3: An exhibition of how the mSPs can be discriminated within SOPs and 
HPs. 



tagged b-jet while the other signature has no tagged b-jets. In this case one finds a 
significant separation of the CPs and HPs from SUPs and SOPs. The lower-left panel 
gives a plot where one signature has two tagged b-jets and the other signature has 
only one tagged b-jet. One again finds that the CPs and HPs are well separated from 
the SOPs and the SUPs for much the same reason as in upper-right panel. Finally, a 
plot is given in the lower-right panel where one signature is the average missing 
while the other signature involves events with no tagged b-jets. Again in this plot the 
CPs (which include mSP4) and HPs are well separated from the SOPs and SUPs. 

The analysis of Fig. fl8.ip exhibits that for some cases, e.g., for the patterns CP and 
HP in the upper right hand corner of Fig. (18. II) . the separation between the SUGRA 
prediction and the Standard Model background is strikingly clear, allowing for the 
identification not only of new physics but also of the nature of the pattern that leads 
to such a signature. 

We discuss now the possibility of discriminating sub-patterns within a given pat- 
tern class. An analysis illustrating this possibility is given in Fig. (18.21) . Here the top 
two panels illustrate how the sub-patterns mSPl, mSP2, mSP4 within the chargino 
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class (CP) are distinguishable with appropriate choice of the signatures. A similar 
analysis regarding the discrimination for the sub-patterns in the stau class (SUP) is 
given in the two bottom panels. The left panel in Fig. (18. 3p gives an analysis of how 
one may discriminate the stop sub-patterns mSPll, mSP12, mSP13 in the stop class 
(SOP), and finally the right panel shows the plots that allows one to discriminate 
the Higgs patterns mSP14 and mSP16 from each other. There are a variety of other 
plots which allow one to discriminate among patterns. With 40 counting signatures 
one can have 780 such plots and it is not possible to display all of them. A global 
analysis where the signatures are simultaneously considered for a large collection of 
mSPs and NUSPs. 

As mentioned in the above analysis we have included models which can produce at 
least 500 SUSY events with 10 fb^^ which is lower than our estimated discovery limits 
for total SUSY events which are about 2200 in this case. The reason for inclusion 
of points below the discovery limit in the total SUSY events is that some of them 
can be detected in other channels such as in the trileptonic channel while others will 
be detectable as the luminosity goes higher. We note in passing that reduction of 
admissible points makes separation of patterns easier. 

8.2 Sparticle Signatures including Nonuniversali- 
ties 

In this section, we give an analysis including nonuniversalities in three different sec- 
tors: NUH, NU3, and NUG. In our analysis we simulate various models with the same 
constant number of events N which we take, as an example to be N= 10^. To dis- 
criminate among the patterns in the signature space, we introduce another set of post 
trigger cuts, which we denote as 'b jet cuts', in addition to the default post trigger 
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Figure 8.4: Discrimination among mSPs within both mSUGRA and NUSUGRA mod- 
els. 

cuts specified in previous chapter. The criteria in the b-jet cuts are the same as the 
default post trigger cuts, except that we change the condition 'at least two hadronic 
jets in the event' to 'specifically at least one b-tagged jet in the event'. We exhibit 
our analysis utilizing both the default cuts and the b jet cuts in Fig.f l8.4p . Two mSPs 
are presented in each figure in different signature spaces to show the separation for 
each case. Signals are simulated with constant number of events in PGS4 for each 
pattern. One can see that even with inclusion of a variety of soft breaking scenarios, 
some mSPs still have very distinct signatures in some specific channels. 

Thus in the top-left panel of Fig. (18. 4p we give a plot of mSP7 (SUP) and mSPll 
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(SOP) in the signature space IL/N (b jet cuts) vs OLlb/N, where OLlb/N is obtained 
with the defauh post trigger cuts. Here we find that these two model types are clearly 
distinguishable as highlighted by shaded and unshaded regions. A similar analysis 
with signatures consisting of ILlb/N (b jet cuts) vs OLlb/N for mSP4 (CP) and mSP7 
(SUP) is given in the top-right panel. The lower-left panel gives an analysis of mSP4 
(CP) and mSP5 (SUP) also in the signature space consisting of ILlb/N (b jet cuts) 
vs OLlb/N. Finally, in the lower-right panel we give an analysis of mSP3 (CP) and 
mSPll (SOP) in the signature plane e"'"/N vs Ib/N. These analyses illustrate that the 
patterns and often even the sub-patterns can be discriminated with the appropriate 
choice of signatures for a general class of SUGRA models including nonuniversalities . 



8.3 The Trileptonic Signal as a Pattern Discrimi- 
nant 
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Figure 8.5: A plot of the number of trilepton events versus the light chargino mass 
for three patterns, one from each class, CP, SUP and HP. The SUP pattern gives the 
largest trileptonic signal followed by the HP and CP patterns. 
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Figure 8.6: The number of tri-lepton events versus the sparticle mass sphttings. 



The trileptonic signal is an important signal for the discovery of supersymmetry. 
For on-shell decays the trileptonic signal was discussed in the early days in flEl 1125] 
and for off-shell decays in |126j . (For a recent application see |118] ). Here we discuss 
the trileptonic signal in the context of discrimination of hierarchical patterns. In 
Fig. (18. 5p we exhibit the dependency of the trilepton signal on the chargino mass. It 
is seen that mSP5 gives the largest number of events in this channel while the CP 
pattern (mSPl) and the HP pattern (mSP14) can also produce a large number of 
trilepton events above the discovery limit, while the chargino mass reach is extended 
for the mSP5 as opposed to the mSPl and mSP14. The above observations hold for 
some of the other SUP patterns as well. Thus the trileptonic signal is strong enough 
to be probed up to chargino masses of about 500 GeV in the SUP pattern. Another 
interesting display of the trileptonic signal is when this signal is plotted against some 
relevant mass splittings. The left panel of Fig. (18.61) shows clear separations for 
hierarchical mass patterns in the number of trilepton events produced with 10 fb~^ 
as a function of the NLSP and the LSP mass splitting for the chargino (CP) pattern 
mSPl and Stau (SUP) mSP5. The plot on the right shows a similar effect for the 
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case where the mass sphtting is taken to be the difference of the CP odd Higgs boson 
mass and the LSP for both the Higgs pattern mSP14 and the stau pattern mSP5. 
The Standard Model background is highly suppressed in this channel. Thus the left- 
panel of Fig. (18. 6p gives an analysis for the trileptonic signal for two patterns: the 
Chargino pattern mSPl and the Stau pattern mSP5 plotted against the NLSP-LSP 
mass splitting with 10 fb~^ of data. 

The analysis of the left-panel of Fig. (18.61) shows that the SUP pattern presents an 
excellent opportunity for discovering SUSY through the 3 lepton mode. The analysis 
also shows a clear separation among mass patterns and further a majority of the 
model points stand above the discovery limit which in this channel is ~ 15 events 
under the post trigger level cuts discussed previously. The right-panel of Fig. (18.61) 
gives an analysis of the trileptonic signal vs the mass splitting of the CP odd Higgs 
and the lightest neutralino LSP for patterns mSP5 and mSP14. Again, we see a 
clear separation of model points. We note that CP odd Higgs can sometimes be even 
lighter than the LSP, and thus the quantity AM = M4 — M^o plotted on the x-axis 
can sometimes become negative. 
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Chapter 9 

Kinematic Signatures of Sparticles 



In this chapter, we discuss various kinematical variables that are useful for pattern 
discrimination. Usually, the kinematic distributions require much higher luminosity 
than the counting signatures, so it is unhkely that SUSY is first discovered with kine- 
matical distributions. But the kinematic distributions carry more precise information 
which is essential for determining the detailed structure of SUSY models. 

Typically, there are two types of kinematic distributions, Pt distributions and 
invariant mass distributions. The Pt distributions we investigate include missing Pt 
distribution and the effective mass distributions. The invariant mass distribution 
analyzed here is the opposite sign same fiavor dilepton invariant mass distribution. 



The kinematical signatures are important for pattern discrimination in addition to 
the event counting signatures discussed previously. We illustrate this using the kine- 
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9.1 Transverse Momentum Distributions 
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matical variables consisting of missing Pt and the effective mass (see Table (17.ip for 
their definitions) and an illustration is given in Fig. (l9.ip . Specifically the analysis of 

1200 1 • • • - 



LHC@10fb"^ 



mSM/500 
mSPI (CP) 
♦ mSP5 (SUP) 
mSP11(S0P) 
+ mSP14(HP) 




Trigger Level Cuts Imposed 



400 600 
Missing P.^. (GeV) 



1000 



Figure 9.1: An exhibition of the missing Pt distributions for 4 different mSUGRA 
models with each corresponding to one class of mSPs, and for the Standard Model. 
Only trigger level cuts are employed here. 



Fig. fl9.ip uses four mSUGRA points one each in the patterns CP, SUP, SOP and HP. 
We exhibit the mSUGRA points used here in the order (mo, mi/2, Aq, tan/?, sign/i): 



CP 


Point 


(3206.9, 


285.3, 


-1319.8, 


9.7, 


+1) 


SUP 


Point 


(92.6, 


462.1, 


352.2, 


4.5, 


+ 1) 


SOP 


Point 


(2296.9, 


625.0, 


-5254.9, 


13.6, 


+1) 


HP 


Point 


(756.8, 


387.0, 


1144.9, 


56.5, 


+1) 



(9.1) 



In the missing P^ distribution the Standard Model tends to produce events with a 
lower missing Pt relative to the mSUGRA case which generates events at relatively 
higher missing Pt- Further, there is a large variation between different mSUGRA 
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Figure 9.2: An exhibition of the effective mass distributions for the same mSUGRA 
models as shown in Fig. (19.11) . Only trigger level cuts are employed here. 

models, as can be seen in Fig. fl9.ip . Thus, for example mSP5 (a stau pattern) and 
mSP14 (a Higgs pattern) have peaks at larger values of missing Pt relative to mSPl 
(a chargino pattern) and mSPll (a stop pattern). Additionally, the shapes of the 
distributions are also different. 

The analysis of effective mass distribution in Fig. (19.21) is carried out with the 
same mSUGRA model points as in Fig. (19. ip . And it is found that in the effective 
mass distribution, the Standard Model tends to produce events with a lower effective 
mass relative to the mSUGRA models, and the variation between mSUGRA models 
remain similar to the case as in Fig. (19.11) . 

The analysis of Fig. (19.11) and Fig. (19.21) shows that the distributions for the 
CP, HP, SOP and SUP are substantially different. It is interesting to note that 
in the missing Pt distribution, the HP and SUP model points have a relatively flat 
distribution compared to the CP and SOP model points. The missing Pt distribution 
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Figure 9.3: The effective mass distributions for 4 different mSUGRA models with 
each corresponding to one class of mSPs, and for the Standard Model. Post trigger 
level cuts are imposed here. The bin size used here is 25 GeV. 



and the effective mass distribution are useful when designing post trigger level cuts to 
optimize the signal over the background. For instance, one can take a 1 TeV effective 
mass cut to analyze the SUP and HP signals shown in Fig. (19. 21) . but this method will 
not work well when it comes to the CP and SOP points since most of their events 
have a rather small effective mass. To illustrate that different models have different 
effective mass distributions, and consequently different effective mass cuts are needed 
for different patterns, an analysis is given in Fig. fl9.3p for the same set of points in 
Fig. fl9.2p with post trigger level cuts imposed. 



9.2 Invariant Mass Distributions 

We also investigate the invariant mass distribution for the opposite sign same flavor 
(OSSF) di-leptons (e+e~, in Fig. (19.41) . We applied the default post trigger 
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cuts as discussed previously to suppress the SM background. As a comparison the 
dominant Standard Model tt background is also exhibited. We have cross checked 
our work with the CMS Note [HE], and found good agreement regarding the SUSY 
signals and the Standard Model background. It is seen that the two mSP points, 



mSP4 Point (1674.9, 137.6, 1986.5, 18.6, +1) 
mSP5 Point (84.4, 429.3, -263, 3.4, +1) 



(9.2) 



plotted in Fig. fl9.4p are clearly distinguishable from each other in the distribution. 
An analysis of invariant mass distribution is given in the Appendix. Here, we apply 

100 1 < < — 
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Figure 9.4: A plot of the opposite sign same flavor (OSSF) di-lepton invariant mass 
distribution at LHC with 10 fb~^ with the default post trigger cuts imposed for two 
different mSP points on top of the SM ti background. 



the general formula given in the Appendix to one specific interesting SUSY decay 
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chains x§ ^ ^ 



For mSP5 model point plotted here, the relevant branching ratios are 



(9.3) 



BR{x^ -^i + £)^ 23.5% and BR{i ^ xl + ^) ^ 100% 



(9.4) 



where i are ei and fii, and ^ are electron and muon. The relevant sparticle masses are 
Mj = 300.9 GeV, M^o = 327 A GeV, and . 
value of the invariant mass of dilepton is 



M^= 300.9 GeV, M^o = 327 A GeV, and M^o = 181.9 GeV. Therefore, the maximum 



M, 



327.4W1 - 



300.92 
327.42 



1 - 



I8I.92 
300.92 



102.8 GeV 



(9.5) 



which is consistent with the result of Fig. (19. 4p . 

For the mSP4 model point plotted here, since the sfermion masses are quite large, 
the SUSY production is dominated by the Ino-production, especially 5^, X2, and X\. 
The decay mode which is responsible for most of the production of the dilepton events 
is 



5.66%. 



(9.6) 



For this model point, the lightest neutralino masses are M^a^ = 54.6 GeV and M^o = 
107 A GeV. Here the decay process is realized through the off-shell slepton decay, 
since the on-shell slepton masses are above TeV. When the off-shell mass of the 



corresponding slepton is moffsheii = ^M^oM^o, the invariant mass edge value achieves 
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its maximum value 



(M^7")^^" = M^o - M^o = 107.4 - 54.6 = 52.8 GeV (9.7) 
which is again consistent with the result of Fig. (19.41) . The off-shell mass of the 




Figure 9.5: Decay process via off-shell slepton. 

slepton m = y^M~oiW~o = V107.4 * 54.6 = 76.6 GeV gives rise to the edge of the 
OSSF dilepton invariant mass distribution. When the off-shell slepton has an off- 
shell mass roughly the half way between the Xi ^"^^ both leptons are likely to 
gain enough transverse momentum in order to pass the detector cuts on lepton Pt- 

The mSP4 model point shown here has recently been investigated |127] in the 
context of helicity amplitudes as a discovery mechanism for supersymmetry. 
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Chapter 10 

Global Analysis of Sparticle 
Patterns: Fuzzy Signature Vectors 



In this chapter, we discuss how one may distinguish sparticle patterns utihzing all 
signature channels that are available. We have given specific examples of how patterns 
can be differentiated from each other. In the previous analysis we used only a few of 
the 40 signatures exhibited in Table (17.11) . Here we want to examined all of them. 

Thus for each parameter point we have analyzed 40 signatures. We now define 
correlations among these signatures. Thus consider an ordered set where the sig- 
natures are labeled 81,82, ..,8^0 and let the number of events in each signature be 
Ni, N2, .., N4Q. Define a signature vector for a given point Xa (a = l,2,..,p) in the 
parameter space 



where = Nf /N and N is the total number of SUSY events. As the parameter point 
Xa varies over the allowed range within a given pattern it generates a signature vector 
where the elements trace out a given range. Thus for a pattern X one generates a 




(10.1) 
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fuzzy pattern vector A^^ so that 

Ae^ = (Aef,Ae2'',..,Aefo), (10.2) 

where A^^ is the range traced out by the element as the parameter point Xa 
moves in the allowed parameter space of the pattern X. What makes the vector A^^ 
fuzzy is that its elements are not single numbers but a set which cover a range. We 
define now the inner product of two such fuzzy pattern vectors so that 

Cxy = (Ae^|AH = 0(l) (10.3) 

where the inner product is if the element A^f and A^f overlap for alH (i = 1, .,40), 
and 1 if at least one of the elements of pattern X, A^^ does not overlap with A^J , 
the element for pattern Y. Therefore, if for two patterns X and Y one finds there is 
no overlap at least for one signature component A^j, then these two patterns can be 
distinguished in this specific signature and one obtains Cxy — 1- Otherwise Cxv — 
which means that all components of A^^ and A^^ have an overlap and cannot be 
distinguished under this criteria. 

We can generalize the above procedure for the signatures 

Cm = ^, (^,J = 1,...,40). (10.4) 

Repeating the previous analysis, one can construct another fuzzy signature vector for 
pattern X as 

AC^ = (AC^„..,AC5,..,AC3^,_,o) (10-5) 
where the elements have a range corresponding to the range spanned by the soft 
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parameters as they move over the parameter space specific to the pattern. Further, 
the definition of the inner product Eq. (110.31) still holds for this new fuzzy signature 
vector. We have carried out a full signature analysis of such comparisons, using 40 
different signatures, and their combinations as defined in Eq. (110.41) and Eq. (110.51) . 





m5 


ml 


m3 


m7 


mil 


m6 


ml2 


ml3 


Nl 


m4 


mis 


m5 














1 





1 


1 





1 





ml 




















1 








1 





m3 














1 





1 


1 











m7 














1 





1 


1 





1 





mil 


1 





1 


1 








1 


1 


1 





1 


m6 




















1 














ml2 


1 


1 


1 


1 


1 


1 








1 


1 


1 


ml3 


1 





1 


1 


1 











1 


1 





Nl 














1 





1 


1 





1 


1 


m4 


1 


1 





1 








1 


1 


1 





1 


ml8 














1 





1 





1 


1 






Table 10.1: A table exhibiting the discrimination of patterns using the criterion of 
Eq. (110.31) where various signatures with both the default post trigger cuts and b jet 
cuts are utilized. If the element of i*^ row and j^^ column is 1, i.e., Cij = 1, one can 
distinguish the z*^ mass pattern from the j^^ one. Here m stands for mSP, and N 
for NUSP. Thus shown here is the discrimination table between the patterns: mSP5, 
mSPl, mSP3, mSP7, mSPll, mSP6, mSP12, mSP13, NUSPl, mSP4, mSPlS. The 
order of the patterns indicates how often these patterns appear in our Monte Carlo 
scan. 



An illustration of the global analysis is given in Tables (110. Ij 110. 2^ 110.31) . We 
carry out this analysis with a large collection of SUGRA model points which belong 
to 22 different hierarchical mass patterns. The complete set of the LHC signatures 
are obtained with the default post trigger cuts as well as the b jet cuts as specified 
in chapter ([8]). We have roughly divided the 22 sparticle patterns into two equal size 
sets, and classified the more probable patterns into the first set: mSP5, mSPl, mSP3, 
mSP7, mSPll, mSP6, mSP12, mSP13, NUSPl, mSP4, mSP18; and the less probable 
patterns into the second set: NUSP13, mSP20, mSPlO, mSP17, NUSP3, mSP19, 
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NUSP5, NUSP8, NUSPIO, NUSP4, NUSP9. The global analysis within these two 
sets is exhibited in Tables (110.11 110.21) . and the analysis between these two sets is 
exhibited in Tables (110. 3p . Altogether, Tables (110.11 110^21 110. 3p show whether or not 
one can distinguish any pair of patterns chosen from the 22 different patterns utilizing 



the signatures investigated here. 
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Table 10.2: A discrimination table between the patterns: NUSP13, mSP20, mSPlO, 
mSP17, NUSP3, mSP19, NUSP5, NUSP8, NUSPIO, NUSP4, NUSP9. 



The analysis shows that it is possible to often distinguish patterns using the cri- 
terion of Eq. (ll0.3p . We note that the analyses exhibited in Fig. (18.41) are the special 
cases of the results in Tables (110. H 110.21 110.5]) . For instance, the clear separation be- 
tween mSP7 and mSPll in the signature space shown in the top-left panel of Fig. (18.41) 
gives the elements C45 = C54 = 1 of Table (110.11) . As indicated in Table (110.21) all 
the patterns analyzed here can be discriminated from each other. This is not really 
surprising, because the probability of finding the sparticle patterns shown in Table 
(110. 2p are not big, and each pattern here is not analyzed with enough model points 
in order to gain sufficient statistics. As emphasized already the analysis of in Tables 
(110. H 110.21 110. 3p is for illustrative purposes as we used a random sample of 22 pat- 
terns out of 37. Inclusion of each additional mass pattern brings in a significant set 
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of model points which need to be simulated, and here one is limited by computing 
power. The full analysis including all the patterns can be implemented along similar 
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Table 10.3: A discrimination table between the two sets of patterns exhibited in the 
previous Tables (fTOTl [To:2|) . 



lines with the necessary computing power. Finally we note that the analysis in Tables 
(110. H 110. 2[ 110. 3p is done without statistical uncertainties. Inclusion of uncertainties 
in pattern analysis would certainly be worthwhile in a future work. 
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Chapter 11 

Signature Degeneracies and 
Resolution of Soft Parameters 



11.1 Lifting Signature Degeneracies 

It may happen that two distinct points in the soft parameter space may lead to the 
same set of signatures for a given integrated luminosity within some predefined notion 
of indistinguishability. Thus consider two parameter points A and B and define the 
'pulls' in each of their signatures by 

p _ \nt - nf\ 

(TAB = \j{8ntf + {5nff + {Snl^f. (11.1) 

Here 5nf ~ sff^ is the uncertainty in the signature events nf, and we estimate the 
SM uncertainty as 5nf^ ~ ^Jy{5nf + 5nf)/2. Here the parameter y parameterizes 
the effect of the SM events, and for the analysis in this section, we take y = 1. In 
other words, if the pull in each of the signatures is less than 5, then the two SUGRA 



85 



parameter space points are essentially indistinguishable in the signature space. In such 
a situation one could still distinguish model points either by including more signatures, 
or by an increase in luminosity. Thus, for example, inclusion of the Higgs production 
cross sections, Bg /^^/^~ constraints, as well as the inclusion of neutralino proton 
scattering cross section constraints tend to discriminate among the model parameter 
points as shown in Ref. PB]. Here we point out that in some cases increasing the 
luminosity can allow one to lift the degeneracies enhancing a subset of signatures in 
one case relative to the other. For illustration we consider the following two sets of 
points in the pattern mSP5 in the mSUGRA parameter space in the following order 
(mo, mi/2, Aq, tan/3, sign/x). 



Point A (192.6,771.3, 1791.1, 8.8, +1), 
Point B (163.0, 761.3, -775.8, 4.7, +1); 



;il.2) 



(11.3) 



Point A' (159.3, 732.3, -783.1, 5.6, +1), 
Point B' (163.5, 753.3, -918.2, 3.3, +1). 

In Table flll.ip we compare the pulls for the pairs of points {A, B) and [A', B') at an 
integrated luminosity of 10 fb~^ and 500 fb~^. For points A and B, one finds that 
the pulls are all less than 2 for an integrated luminosity of 10 fb^^. However, for 
an integrated luminosity of 500 fb~^, the pulls for signatures (5,7,8,34,39) increase 
significantly and the pull for signature number 34 is in excess of 5 allowing one to 
discriminate between the two parameter points A and B. A very similar analysis 
is carried out for parameter points A' and B'. Here one finds that the signature 
(0, 1, 2, 3, 6, 12, 14, 32, 33, 36) receive a big boost as we go from 10 fb^^ to 500 fb~^, 
and the signatures (0, 1, 2, 6, 36) give pulls greater than 5, with the largest pulls being 
in excess of 12, allowing one to discriminate between the parameter points A' and 
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B'. We note the analysis ignores systematic errors and also does not consider an 
ensemble of simulations. Nonetheless it does illustrate the effects of moving from a 
low to a high LHC luminosity allowing one to discriminate some model pairs, which 
appear degenerate in the signature space at one luminosity, but can become distinct 
from each other at a larger luminosity. 



11.2 Resolving Soft Parameters using LHC Data 
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Figure 11.1: An analysis showing the resolutions in mo and mi/2 that can be reached 
with 1000 fb~^ of integrated luminosity under the REWSB constraints. The two left 
panels give the number of SUSY events vs mo (top left panel) and vs mi/2 (lower left 
panel) for 1000 fb~^ of integrated luminosity. The right panels give the resolutions in 
mo (top right panel) and in mi/2 (lower right panel) using the left panels. 
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We discuss now the issue of how well we can resolve the points in the parameter 
space Xq, (a = 1, for a given luminosity. Consider Eg. (11 1.1 1) and set 5N = \fN , 
and parameterize the standard model uncertainty by 5N^^'^ = ^Jy5N . Next we set 
the criterion for the resolution of two adjacent points in the SUGRA parameter space 
separated by Ax^ so that the separation in the signature space satisfies 



Since = (Ts\isy{,Xa)Cu:ic^ where (Tsusy is the cross section for the production of spar- 
ticles, and £lhc is the LHC integrated luminosity, the resolution achievable in the 
vicinity of SUGRA parameter point Xa at that luminosity is given by 



In Fig. (Ill.ip we give an illustration of the above when itlq varies between 500 GeV and 
2000 GeV while mi/2 = 500 GeV, Aq = 0, tan/5 = 30, and > 0. From Fig. lITLTl) one 
finds that the resolution in niQ strongly depends on the point in the parameter space 
and on the luminosity. Quite interestingly a resolution as small as a few GeV can be 
achieved for mo in the range 500-1000 GeV with 1000 fb~^ of integrated luminosity. A 
similar analysis varying mi/2 in the range 500~900 GeV for the case when tjiq = 500 
GeV, Aq = 0,tan/5 = 30 and /i > 0, shows that a resolution in mi/2 as low as 1 GeV 
can be achieved with 1000 fb~^ of integrated luminosity. 
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Table 11.1: An exhibition of lifting the degeneracy of two points in the mSUGRA 
parameter space using luminosity. Two pairs of points (^4, B) and (A', B') are 
indistinguishable under the 2 sigma criteria at 10 fb~^ luminosity (column 3-8), but 

can be clearly separated when the himinosity increases to 500 fb~^ (column 9-14). 
The Standard Model uncertainty is estimated as 5nf^ = {5nf -\- 5nf)/2. 
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Chapter 12 

Higgs Production at Colliders 



In this chapter we investigate the Higgs cross sections at the Tevatron and at the 
LHC The hghtness of A (and also of H and iJ^) in the Higgs Patterns implies that 
the Higgs production cross sections can be large (for some of the previous analyses 
where light Higgses appear see |128[ I129[ 1104] ) . 

Quite interestingly the recent Tevatron data is beginning to constrain the Higgs 
Patterns (HPs). This is exhibited in the left panel of Fig. (112. 11) where the leading 
order (LO) cross section for the sum of neutral Higgs processes a^rripp) = [(^{pp 
$)BR($ 2t)] (where sum over the neutral $ fields is implied) vs the CP odd 
Higgs mass is plotted for CM energy of y/s = 1.96 TeV at the Tevatron. One finds 
that the predictions of cqtrripp) from the HPs are the largest and lie in a narrow 
band followed by those from the Chargino Pattern mSP2. The recent data from the 
Tevatron is also shown [130] . A comparison of the theory prediction with data shows 
that the HPs are being constrained by experiment. Exhibited in the right panel of 
Fig. (112. II) is (Jq>TT{pp) = [(^{pp $)BR($ 2r)] arising from the HPs (and also 
from other patterns which make a comparable contribution) vs the CP odd Higgs 
mass with the analysis done at CM energy of ^/s = 14 TeV at the LHC. Again it is 
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Figure 12.1: Left panel: Predictions for [a{pp $)BR($ 2r)] in mSUGRA as a 
function of the CP odd Higgs mass tua for the HPs at the Tevatron with CM energy 
of y/s = 1.96 TeV. The hmits from D0 are indicated |130j . Right paneh Predictions 
for [a{pp $)BR(<I> ^ 2r)] in mSUGRA as a function of tua at the LHC with CM 
energy of y/s = 14 TeV for the HPs, the chargino pattern mSPl and the stau pattern 
mSP5. The HPs are seen to give the largest cross sections. 



seen that the predictions of cr$rr(pp) arising from the HPs are the largest and lie in a 
very narrow band and the next largest predictions for cr$^T-(pp) are typically from the 
Chargino Patterns (CPs). The larger cross sections for the HPs enhance the prospects 
of their detection. 

Since the largest Higgs production cross sections at the LHC arise from the Higgs 
Patterns and the Chargino Patterns we exhibit the mass of the light Higgs as a 
function of rriQ for these two patterns in the left panel of Fig. (112.21) . We note that 
many of the Chargino Pattern points in this figure appear to have large mo indicating 
that they originate from the Hyperbolic Branch/Focus Point (HB/FP) region 

We discuss now briefly the Higgs to bb decay at the Tevatron. From the pa- 
rameter space of mSUGRA that enters in Fig.(l) we can compute the quantity 
[{pp $)BR($ — s> bb)]. Experimentally, however, this quantity is difficult to mea- 
sure because there is a large background to the production from qq, gg bb. For 
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Figure 12.2: Left panel: mSPl and HPs are plotted in the mo-mh plane in mSUGRA 
fi > 0. Right panel: Predictions for [cr(pp <I>)BR($ 2r)] in NUSUGRA 
(NUH,NUG,NU3) as a function of CP odd Higgs mass at the LHC showing that 
the HPs extend beyond 600 GeV with non-universalities (to be compared with the 
analysis of Fig. (112. II) under the same naturalness assumptions). 



this reason one focuses on the production [{pp <I>6)BR($ bb)] \131\ . For the 
parameter space of Fig.(l) one gets [{pp $6)BR($ bb)] < 1 pb at (tan/3 = 
55, Ma = 200 GeV). The preliminary CDF data [132] puts limits at 200 GeV, in the 
range (5-20) pb over a 2a band at the tail of the data set. These limits are larger, 
and thus less stringent, than what one gets from $ t^t^ . For the LHC, we find 
[{pp $6)BR($ hb)] ~ 200 pb for the same model point. A more detailed fit 
requires a full treatment which is outside the scope of the present analysis. 

The neutral Higgs production cross section for the NUSUGRA case is given in the 
right panel of Fig. (112.21) . The analysis shows that the Higgs Patterns produce the 
largest cross sections followed by the Chargino Patterns as in mSUGRA case. One 
feature which is now different is that the Higgs Patterns survive significantly beyond 
the CP odd Higgs mass of 600 GeV within our assumed naturalness assumptions. 
Thus nonuniversalities tend to extend the CP odd Higgs beyond what one has in the 
mSUGRA case. 
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Chapter 13 

Bg IJ-~^I^~ Constraints 

In this chapter, we investigate the Bg —>■ constraints within the context of 

the sparticle pattern analysis. The process Bs —>■ fJ'^fi~ is dominated by the neutral 
Higgs exchange [133j . The decay B^, £^£~ {d'=d,s) is governed by the effective 
Hamiltonian [134] 

Hes = -£^VtbV;ACsOs + CpOp + C'sO's + C'pO'p + CloOlo) (13.1) 

where 



Os = mtid'^PRb^)ii£), (13.2) 

Op = m,{d'^PRho){h^£), (13.3) 

O's = ma\d'^PLha)m, (13.4) 

O'p = md^id'^PMih^i), (13.5) 

Oio = KYPLbaKhf^l^i) (13.6) 
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The branching ratio B{B^, 



is then given by 



Ami 



'1- 



\fs? + \fp + 2m,fA\ 



where /j (i=S,P) and /a are defined as follows 



2f B^l y . 
nidi + 



2M 



-C 



10- 



Specifically Cs and Cp have the form 



3 2 



v^m2^cos3/3^^ >^^Mi 



irti tan^ /9 



EE 



;i3.7) 



;i3.8) 
;i3.9) 



13.10) 



13.11) 



\/2m^ cos/3 

When tan (3 becomes large, one finds that the branching ratio is propotional to tan^ (3. 

It is thus reasonable to expect that the Higgs patterns (HPs) will be constrained 
more severely than other patterns by the Bs experiment, since HP points 

usually arise from the high tan (3 region (however, as we noticed already, the nonuni- 
versalities in the Higgs sector (NUH) can also give rise to HPs for moderate values of 
tan [3). 

In Fig. (113. II) we carry out a detailed analysis where the branching ratio Br{Bs 
is plotted against the CP odd Higgs mass m^. The upper left (right) hand 
panel gives the analysis for the case of mSUGRA for /i > (/i < 0) for the Higgs 
Patterns as well as for several other patterns, and the experimental constraints are 
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also shown. One finds that the constraints are very effective for ^ > (but not for 
fi < 0) constraining a part of the parameter space of the HPs and also some models 
within the Chargino and the Stau Patterns are constrained (see upper left and lower 
left panels of Fig. ffTO) ). 

From the analysis of Fig. (113. II) . it is observed that the strict imposition of the 
constraint Br{Bs Ai+yU^) < 1.5 x 10^'' still allows for large tan/? in the mSUGRA 
model. Thus all of the HP model points given in Fig. (113. II) that satisfy this constraint 
for the mSUGRA fi > case correspond to tan/5 in the range of 50 - 55. A similar 
limit on tan j3 is also observed for the nonuniversal models. We remark, however, that 
the HPs are not restricted to large tan j3 in particular for the case of the NUH model. 
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Constraints : n h e WMAP, Br{b^ s y), and mass limits 
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Figure 13.1: Predictions for the branching ratio Bg in various patterns 

in the SUGRA landscape. Upper left panel: predictions are for the patterns for 
/i > in mSUGRA; upper right panel: predictions are for the patterns for < 
in mSUGRA; lower left panel: predictions for the Higgs Patterns alone for /i > in 
mSUGRA; lower right panel: predictions for NUSUGRA models NUH, NUq3, and 
NUG for /i > 0. The experimental limits are: top band 2005 [1351 1136j . and the 
bottom two horizontal lines are preliminary limits from the CDF and D0 data |77j . 
For convenience we draw the limits extending past the observable mass of the CP 
odd Higgs at the Tevatron. 
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Chapter 14 

Direct Detection of Dark Matter 

14.1 mSUGRA 



We discuss now the direct detection of dark matter within the framework of the 
mSUGRA models. In direct detection experiments one measures the cross section of 
the WIMP scattering off the heavy nuclei such as germanium. The neutralino inter- 
acts with quarks in the target nuclei through a Higgs boson exchange, or a squark ex- 
change. The neutralino-nucleus scattering cross-section typically is dominated by the 



scalar part of the neutralino-quark interaction and thus it is the quantity a^p{scalar) 
that is of interest to us. The basic interaction governing the x ~ P scattering is the 
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effective four-fermi interaction given by (see e.g. |137j ) 



ayp(scalar) = 

TT 



E 1- E /n E 

i=uA,s \ i=uA,s I a=c,b,t 



Aff = Xliil^XlYi^Pi + 5PR)g + Cxxmqqq + Dx75X"^g9755 

+Exil5Xmgqq + Fxxmgqn^q. (14.1) 

Tlie X — p cross-section arising from scalar interactions Cxx^T^qQQ is given by 

(14.2) 

Here fir is tlie reduced mass, /f (i=u,d,s quarks) are defined by 

rupff =< p\mgiqiqi\p >, (14.3) 

and C is given by 

C = Cho + Cho + C/, (14.4) 

wliere C/jO,C//o are tlie contributions from the s-channel and exchanges and 
is the contribution from the t-channel sfermion exchange. They are given by [137| 

C>.o{u,d) = -{+)-rJ^'^^4r^Rea. (14.5) 

CMu,d) = (14.6) 

4Mvi/iVf^o sm p(cosp) 

Here (u,d) refer to the quark flavor, a is the Higgs mixing angle, and Cgi, C^j^ etc. 
are as defined in Ref . [137] , and a and p are defined by 

cr = X4q(X2o — tan 6*1^-^10) cos a + X^q{X2q — ta.n9wXlQ) sin a, (14-8) 
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p = —XIq (XgQ — tan Ow^io) ^ + ^30(^^20 ~ ^w^io) cos a, (14.9) 
where X„o are the components of the LSP 

X = x*,B + x;,Ws + x;,H, + xi^h^. (m.io) 

The coefficients /f are associated with some amount of uncertainties |138l I129j 

= 0.020 ±0.004, (14.11) 
fl = 0.026 ±0.005, (14.12) 
fP = 0.118 ±0.062. (14.13) 

In the absence of CP phases, in Fig. (114. II) [23j we give an analysis of the scalar 
neutralino-proton cross section cr{xiP) as a function of the LSP mass (for a sample 
of Post- WMAP3 analysis of dark matter see fl9[ I139j , and for more recent analysis 
see [231 1140[ I14H [Ti2j ). The upper left panel of Fig. (114. II) gives the scalar a{xiP) for 
the mSUGRA parameter space for p > 0. We note that the Higgs patterns typically 
give the largest dark matter cross sections (see the upper left and lower left panels 
of Fig. (114. 11) ) and are the ffist ones to be constrained by experiment. The second 
largest cross sections arise from the Chargino Patterns which shows an embankment, 
or Wall, with a copious number of points with cross sections in the range 10~^^^'^cm^ 
(see the upper left panel and lower right panel), followed by Stau Patterns (lower left 
panel), with the Stop Patterns producing the smallest cross sections (upper left and 
lower right panels). The upper right panel of Fig. (114. II) gives the scalar cross section 
cr{xiP) for yU < and here one finds that the largest cross sections arise from the 
CPs which also have a Chargino Wall with cross sections in the range 10~^^^'^cm^ 
(upper right panel). The analysis shows that altogether the scalar cross sections he 
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Figure 14.1: Analysis of oixv) for mSUGRA: upper left panel: > case including 
all patterns; upper right panel: fi < allowing all patterns; lower left hand panel: A 
comparison of cr{xp) ^oi HPs and a stau NLSP case which is of type mSP5 for /i > 0; 
lower right panel: a comparison of cr(xp) for the Chargino Pattern mSPl vs the 
Stop Patterns mSPll-mSP13. The analysis shows a Wall consisting of a clustering 
of points in the Chargino Patterns mSPl-mSP4 with a (t{xp) the range io~^^^-^ 
cm^ enhancing the prospects for the observation of dark matter by SuperCDMS [143j . 
ZEPLIN-MAX[lll] or LUX[Tl5] in this region. 
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in an interesting region and would be accessible to dark matter experiments currently 
underway and improved experiments in the future [1461 11471 11481 11491 11501 I145j . 
Indeed the analysis of Fig. (114. ip shows that some of the parameter space of the Higgs 
Patterns is beginning to be constrained by the CDMS and the XenonlO data |150j . 

14.2 Nonuniversalities of Soft Breaking 

As already discussed in previous chapters, it is useful to consider other soft breaking 
scenarios beyond mSUGRA, since the nature of physics at the Planck scale is largely 
unknown. One such possibility is to consider nonuniveralities in the Kahler potential, 
which can give rise to nonuniversal soft breaking consistent with flavor changing 
neutral current constraints. We consider three possibilities which are nonuniversalities 
in (i) the Higgs sector (NUH), (ii) the third generation squark sector (NU3), and (iii) 
the gaugino sector (NUG) (for a sample of previous work on dark matter analyses with 
nonuniversalities see |151j ). We parametrize these nonuniversalities as in Eq. (16. 2p . 
In each case we carry out a Monte Carlo scan of 1 x 10^ models. The above covers a 
very wide array of models. 

The analysis of the direct detection of dark matter in NUSUGRA are presented in 
Fig. (114.21) . As in the mSUGRA case one finds that the largest dark matter cross sec- 
tions still arise from the Higgs Patterns followed by the Chargino Patterns within 
the three types of nonuniversality models considered: NUH (upper left panel of 
Fig.dH^D), NU3 (upper right panel of Fig.([Ti3)), NUG (lower panel of Fig.([Ti2D). 
Again the analysis within NUSUGRA shows the phenomenon of the Chargino Wall, 
i.e., the existence of a copious number of Chargino Patterns (specifically mSPl) in all 
cases with cross sections in the range 10~'^^^'^cm^. Most of the parameter points along 
the Chargino Wall lie on the Hyperbolic Branch/Focus Point (HB/FP) region 
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where the Higgsino components of the LSP are substantial (for a review see 
Thus this Chargino Wall presents an encouraging region of the parameter space where 
the dark matter may become observable in improved experiments. 

It is seen that Higgs Patterns (HPs) arising in a wide range of models: in mSUGRA, 
and in NUSUGRA models are typically seen to lead to large Higgs production cross 
sections at the Tevatron and at the LHC. It is also seen that the HPs lead typically 
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Figure 14.2: Analysis of the scalar cross section oixv) for NUSUGRA models: NUH 
(upper left panel), NU3 (upper right panel), NUG (lower panel). As in Fig. (114. 11) 
the Wall consisting of a clustering of points in the Chargino Patterns mSPl-mSP4 
persists up to an LSP mass of about 900 GeV with a oixv) in the range io~^^^-^ 
cm^ enhancing the prospects for the observation of dark matter by SuperCDMS and 
ZEPLIN-MAX in this region. 



102 



to the largest neutralino-proton cross sections and would either be the first to be 
observed or the first to be constrained by dark matter experiment. The analysis pre- 
sented here shows the existence of a Chargino Wall consisting of a copious number of 
parameter points in the Chargino Patterns where the NLSP is a chargino which give a 
(t{xiP) at the level of 10~^^^'^cm^ in all models considered for the LSP mass extending 
up to 900 GeV in many cases. These results heighten the possibility for the observa- 
tion of dark matter in improved dark matter experiments such as SuperCDMS jl43j . 
ZEPLIN-MAX[T11], and LUX[Tl5] which are expected to reach a sensitivity of 10""^^ 
cm^ or more. Finally, we note that several of the patterns are well separated in the 
(j(xiP)- LSP mass plots, providing important signatures along with the signatures 
from colliders for mapping out the sparticle parameter space. 
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Chapter 15 



Conclusions 



The minimal supersymmetric Standard Model has 32 sparticle masses. Since the soft 
breaking sector MSSM is arbitrary, one is led to a landscape of as many as 10^^ or 
more possibilities for the sparticle mass hierarchies. The number of possibilities is 
drastically reduced in well motivated models such as supergravity models, and one 
expects similar reductions to occur also in gauge and anomaly mediated models, and 
in string and branc models. We have analyzed the mass hierarchies for the first four 
lightest sparticle (aside from the lightest Higgs boson) for supergravity models. Specif- 
ically, we analyzed the mass hierarchies for the mSUGRA model and for supergravity 
models with nonuniversalities in the soft breaking in the Higgs sector, nonuniversal- 
ities in the soft breaking in the third generation sector, and nonuniversalities in the 
soft breaking in the gaugino sector. It is found that in each case only a small number 
of mass hierarchies or patterns survive the rigorous constraints of radiative breaking 
of the electroweak symmetry, relic density constraints on cold dark matter from the 
WMAP data, and other experimental constraints from colliders. These mass hier- 
archies can be conveniently put into different classes labeled by the sparticle which 
is next heavier after the LSP. For the SUGRA models we find six different classes: 
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chargino patterns, stau patterns, stop patterns, Higgs patterns, neutralino patterns, 
and gluino patterns. 

We discussed the techniques for the analysis of the signatures and the technical 
details on simulations of sparticle events. We also discuss the backgrounds to the 
SUSY phenomena arising from the Standard Model processes. Additionally we dis- 
cussed the identification of patterns based on 40 event identification criteria. It is 
found that these criteria allow one to discriminate among most of the patterns. An 
analysis of how one may lift degeneracies in the signature space, and how accurately 
one can determine the soft parameters using the LHC luminosities is also given. In 
addition, we also investigate the Higgs production at the Tevatron and at the LHC, 
and the direct detection of dark matter within the context of the sparticle pattern 
analysis. 

It is hoped that the analyses of the type discussed here would help not only in the 
search for supersymmetry but also allow one to use the signatures to extrapolate back 
to the underlying supersymmetric model using the experimental data when such data 
from the LHC comes in. In the above our analysis was focused on supergravity unified 
models. However, the techniques discussed here have a much wider applicability to 
other models, including models based on gauge and anomaly mediated breaking, as 
well as string and brane based models. 
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Appendix 



Dilepton Invariant Mass 

Here we give some further details of the analysis of the kinematic signatures dis- 
cussed in chapter 1^. Specifically we study the kinematics of the dilepton invariant 
mass of SUSY chain decays. The process we consider here has two successive decays, 



Figure 15.1: Particle A decays into particle B and C, and particle C continues to 
decay into particle D and E. 

A ^ B + C followed hj C ^ D + E as shown in Fig. (115.11) . Let particles B and 
D be the Standard Model particles and the A, C and E be the SUSY particles. We 
will make the approximation that all SM particles are massless (except for the top 
quark) . 

Let us consider the decay process A ^ B + C m the rest frame of A. We will adopt 
the following notation A = rriA, B = itlb, and C = me- Using energy-momentum 
conservation, one can obtain the following relations 



A 



B 




D 



E 



P 



^A^ -{B + Cy^A^ -{B- C)^ 
2A 



(15.1) 



Eb 



Ec 



2A 



(15.2) 



(15.3) 
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where P = \Pb\ = \Pc\- We can do the same calculation for the second decay process 
C ^ D + E in the rest frame of C by the following substitutions A ^ C, B ^ D, 
and C E. However the calculation for these two successive decay process is carried 
out in two different inertial frames. Thus we introduce the following notation (-P^)a 
to indicate that the four vector momentum of particle B is defined in the rest frame 
of particle A. So we rewrite the Eqs. fll5.mi5.HI) ioi A ^ B + C and C ^ D + E, 
and take the approximation that B = D = 0. 

iPB,c)A = ^^, {Es)a = ^^, {Ec)a = ^^; (15.4) 

(^2 _ p2 fi2 _ p2 fi2 I 7712 

{Pd,e)c = ^^^^, iEo)c = ^^, {Ee)c = ^^. (15.5) 

However, in order to reconstruct the invariant mass for the Standard Model parti- 
cles B and D, we have to obtain the energy- momentum four- vector for both particles 
in the one frame. Thus one performs some Lorentz transformations to convert the 
energy-momentum vectors to the same frame, for instance, transforming {Pd)c to 
{Pd)a- To do this, one has to know {Pc)a which has been done in Eq. (115. 4p 

iEc)A = ^^. (Pc)a = ^j^. (15.6) 
And the Lorentz transformations cirG clS follows 



{PE)a = smOiPoh, (15.7) 
{Pd)a = {lc)A[cos9{Pn)c + mA{En)c], (15.8) 
{Ed)a = {lc)A[{ED)c + cose{Pn)c{f3c)A], (15.9) 

where the angle 9 is the angle between momentum {Pd)c and {Pc)a as shown in 
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Fig. (I15.2P : sm9{PD)c and cos9{Pd)c are the transverse and longitudinal components 
of the momentum of particle D in the rest frame of particle C; and the Lorentz 
transformation variables 7 and f3 are defined as 

(7oU = %^, (Pc)a-^. (15.10) 



Figure 15.2: Angle 6 is the angle between momentum {Pd)c and {Pc)a- The z 
direction, or the longitudinal direction in the Eqs. fll5.7m5.9|) is along the direction 
of {Pc)a. 



In the rest frame of particle A, the Lorentz invariant quantity, the invariant mass 
of the Standard Model particles B and D now can be calculated easily so that 



MsD = ^{Eb + En? - (PEr - (PE - PhY (15.11) 

where we have dropped the subscript A in Eq. (115.111) and will do so in the subsequent 
analysis. The Eq. (115. lip can be simplified 



l + cosO 

Mbd = A\II-—^1- ^ V (15.12) 



which has a maximum value when cos 6' = 1 so that 



= A^/l-^^l-g, (15.13) 
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and it has a vanishing minimum when cos 6 = —1. The results discussed here are 
utihzed in the analysis given in chapter Q. 
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Benchmarks 
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50.0 


-1- 


(0,0) 


(0,0) 


(0,0) 


mSP2 


2365 


1395 


3663 


42.2 




(0.0) 


(0.0) 


(0.0) 


mSP2 


13(>5 


r)!)5 


3012 


35.1 


+ 


(0.110.-0.338) 


(0.0) 


(0,0) 


mSP2 


1166 


507 


-954 


59.6 


+ 


(0,0) 


(0.325,0.458) 


(0,0) 


mSP2 


1414 


221 


-551 


54.3 


+ 


(0,0) 


(0,0) 


(0.156,0.968) 


mSP3 


741 


551 





50.0 


-1- 


(0,0) 


(0,0) 


(0,0) 


mSP3 


1585 


1470 


3133 


39.1 




(0,0) 


(0,0) 


(0,0) 


mSP3 


694 


674 


-1564 


27.0 


+ 


(0.922,-0.293) 


(0,0) 


(0,0) 


mSP3 


570 


559 


1042 


41.3 


-1- 


(0,0) 


(-0.482,-0.202) 


(0,0) 


mSP3 


392 


312 


320 


41.3 


+ 


(0,0) 


(0,0) 


(-0.404,0.908) 


mSP4 


1674 


137 


1985 


18.6 


+ 


(0,0) 


(0,0) 


(0,0) 


mSP4 


1824 


127 


-1828 


6.4 




(0,0) 


(0,0) 


(0,0) 


mSP4 


1021 


132 


-638 


6.6 


+ 


(0,0) 


(-0.020,0.963) 


(0,0) 


mSP4 


2181 


127 


-3859 


3.9 


+ 


(0,0) 


(0,0) 


(0.836,-0.248) 


NUSPl 


2738 


1689 


-4243 


42.4 


+ 


(0,0) 


(-0.828,-0.899) 


(0,0) 


NUSPl 


540 


1190 


2516 


13.9 


+ 


(0,0) 


(0,0) 


(-0.408,-0.660) 


NUSP2 


845 


726 


-75 


48.4 


+ 


(0,0) 


(-0.694,-0.400) 


(0,0) 


NUSP3 


396 


1018 


-179 


18.3 


+ 


(0,0) 


(0,0) 


(0.250,-0.452) 


NUSP4 


400 


1558 


2511 


5.9 


-1- 


(0,0) 


(0,0) 


(-0.401,-0.607) 



Tabic 15.1: Benchmarks for the class CP where the chargino Xi is the NLSP in 
mSUGRA and in NUSUGRA models. Benchmarks are computed with mi^^{m},) — 
4.23 GeV, a^{Mz) = .1172, and mt(pole) - 170.9 GeV. 
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Stan Patterns (Si;Ps) 



SUGRA 


mo 


mi/2 




tan /3 




NUH 


NU3 


NUG 


Pattern 










(sign) 




\0q?,,0thR) 


\pM-2 , OM3 ) 


mSP5 


111 


531 





5.0 


-1- 


(0,0) 


(0,0) 


(0,0) 


mSP5 


162 


569 


1012 


15.8 




(0.0) 


(0.0) 


(0.0) 


iiiSP5 


191 


5 15 


-722 


17.2 


+ 


(-0.3 10.-0.332) 


(0,0) 


((),()) 


mSP5 


114 


440 


-50 


15.2 


+ 


(0,0) 


(-0.204,-0.846) 


(0,0) 




/ 


o^o 


OUi 




-r 




li\ ^^ 


l^u.zo^,-u.uoy ) 


mSP6 


245 


370 


945 


31.0 


-1- 


(0,0) 


(0,0) 


(0,0) 


mSP6 


1452 


1651 


2821 


38.5 




(0,0) 


(0,0) 


(0,0) 


mSP6 


356 


545 


927 


31.7 


+ 


(0.667,0.055) 


(0,0) 


(0,0) 


mSP6 


442 


463 


1150 


11.0 


+ 


(0.0) 


(-0.187.-0.546) 


(0.0) 




■ )Uo 






.>!>. I) 








/ 11 ".i^".: ri 

l^-U..^n.^.U. iU-j f 


mSP7 


75 


201 


230 


14.0 


+ 


(0,0) 


(0,0) 


(0,0) 


mox^ / 


iOL 




yoo 


oD.o 




/n tw 
l^u,u ] 


(^U,U ) 


I^U,U ) 


mSP7 


428 


671 


484 


43.8 


+ 


(-0.392,-0.808) 


(0,0) 


(0,0) 


mSP7 


226 


426 


944 


27.1 


+ 


(0.0) 


(0.176.-0.430) 


(0.0) 


rnSP7 


li:i 


125 


2()G 


23. 1 


+ 


(0.0) 


(0.0) 


(0.718.0.100) 


mSP8 


1880 


877 


4075 


54.8 


+ 


(0,0) 


(0,0) 


(0,0) 


mSP8 


994 


1073 


3761 


38.1 




(0,0) 


(0,0) 


(0,0) 


mSP8 


602 


684 


805 


49.6 


+ 


(0.490,0.326) 


(0,0) 


(0,0) 


mSPS 


470 


624 


-88 


55.4 


+ 


(0,0) 


(-0.531,-0.075) 


(0,0) 


mSPS 


525 


450 


642 


56.4 


+ 


(0,0) 


(0,0) 


(0.623,0.246) 


mSP9 


667 


1154 


-125 


51.0 


+ 


(0,0) 


(0,0) 


(0,0) 


mSP9 


560 


1156 


-1092 


39.5 


_ 


(0,0) 


(0,0) 


(0,0) 


mSP9 


362 


602 


268 


37.0 


+ 


(0.969,-0.232) 


(0,0) 


(0,0) 


mSP9 


496 


731 


679 


49.3 


+ 


(0,0) 


(-0.241,-0.452) 


(0,0) 


mSP9 


485 


478 


-128 


52.8 


+ 


(0,0) 


(0,0) 


(0.971,0.653) 


mSPlO 


336 


772 


-3074 


10.8 


+ 


(0,0) 


(0,0) 


(0,0) 


mSPlO 


738 


1150 


-4893 


15.5 


+ 


(0,0) 


(0.802,0.343) 


(0,0) 


mSP17 


908 


754 


5123 


25.4 




(0,0) 


(0,0) 


(0,0) 


mSPlS 


344 


686 


-2718 


13.8 




(0,0) 


(0,0) 


(0,0) 


mSPlS 


322 


806 


-3069 


9.3 


+ 


(0.526,-0.707) 


(0,0) 


(0,0) 


mSPlS 


60 


290 


-339 


5.2 


+ 


(0,0) 


(0,0) 


(0.967,-0.074) 


mSPlQ 


1530 


1875 


13081 


16.3 




(0,0) 


(0,0) 


(0,0) 


mSP19 


1828 


1326 


-5102 


32.3 


+ 


(0.592,-0.213) 


(0,0) 


(0,0) 


mSP19 


782 


637 


2688 


37.9 


+ 


(0,0) 


(0,0) 


(0.451,-0.551) 


NUSP5 


649 


955 


-1984 


33.5 


+ 


(0,0) 


(-0.763,0.701) 


(0,0) 


NUSP6 


1360 


1736 


-2871 


46.1 


+ 


(0,0) 


(-0.466,0.694) 


(0,0) 


NUSP7 


1481 


1531 


-3169 


42.2 


+ 


(0,0) 


(0,0) 


(0.117,-0.463) 


NUSP8 


670 


1788 


371 


57.9 


+ 


(0,0) 


(0,0) 


(-0.223,0.931) 


NUSP9 


46 


1938 


-48 


13.0 


+ 


(0,0) 


(0,0) 


(-0.412,-0.650) 



Table 15.2: Benchmarks for the class SUP where the stau ti is the NLSP in mSUGRA 
and in NUSUGRA. 
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Stop Patterns (SOPs) 



SUGRA 

Pattern 


mo 

(GeV) 


^1/2 

(GeV) 


(GcV) 


tan /3 


M 

(sign) 


NUH 


NU3 


NUG 

(^Af2 ' ^Ai^s ) 


v>iG;t3i 1 
iTlox^ ± ± 


OIL 




-4o00 


inn 
lU.U 


+ 




(r\ n^ 
(^U,U ) 


{n n^ 


mor^ ± ± 


iOOo 




10 1^ 


K n 
o.y 






(r\ n^ 
(^U,U ) 


{n n^ 


mSPll 


1391 


1089 


8192 


14.9 




(f) 470 632^ 


(0,0) 


fO 0) 


mSPll 


2204 


933 


-1144 


35.6 


-1- 


(0,0) 


(0.642,-0.400) 


(0,0) 


mSPll 


1406 


1471 


-2078 


8.3 


+ 


(0,0) 


(0,0) 


(-0.130,-0.690) 


mSP12 


1371 


1671 


-6855 


10.0 


+ 


(0,0) 


(0,0) 


(0,0) 


mSP12 


1054 


1372 


-5754 


13.7 




(0,0) 


(0,0) 


(0,0) 


mSP12 


915 


927 


-3993 


20.7 


+ 


(0.078,0.833) 


(0,0) 


(0,0) 


mSP12 


826 


1016 


-3926 


12.8 


+ 


(0,0) 


(-0.630,-0.490) 


(0,0) 


mSP12 


1706 


1287 


-4436 


29.7 


+ 


(0,0) 


(0,0) 


(0.416,-0.260) 


mSP13 


524 


800 


-3315 


15.0 


+ 


(0,0) 


(0,0) 


(0,0) 


mSP13 


765 


1192 


-4924 


12.0 




(0,0) 


(0,0) 


(0,0) 


mSP13 


1055 


1601 


-6365 


13.6 


+ 


(0.277,-0.820) 


(0,0) 


(0,0) 


mSP13 


1073 


1664 


-6528 


11.6 


+ 


(0,0) 


(0.728,0.060) 


(0,0) 


mSP13 


540 


774 


-2432 


5.3 


+ 


(0,0) 


(0,0) 


(0.705,-0.201) 


mSP20 


1754 


840 


7385 


13.3 




(0,0) 


(0,0) 


(0,0) 


mSP21 


792 


845 


6404 


12.6 




(0,0) 


(0,0) 


(0,0) 


NUSPIO 


718 


467 


1657 


19.0 


+ 


(0,0) 


(0,0) 


(0.023,-0.810) 



Table 15.3: Benchmarks for the class SOP where the stop ti is the NLSP in mSUGRA 
and in NUSUGRA models. 



Higgs Patterns (HPs) 



SUGRA 
Pattern 


mo 
(GcV) 


mi/2 
(GeV) 


Ao 
(GeV) 


tan P 

{■Vu/Vd) 


(sign) 


NUH 


NU3 

(5<33> ^tbfl) 


NUG 
{^M2 1 ^Ma ) 


mSP14 


1040 


560 


450 


53.5 


+ 


(0,0) 


(0,0) 


(0,0) 


mSP14 


760 


515 


2250 


31.0 


+ 


(0.255,-0.500) 


(0,0) 


(0,0) 


mSP14 


740 


620 


840 


53.1 


+ 


(0,0) 


(-0.530,-0.249) 


(0,0) 


mSP14 


1205 


331 


-710 


55.0 


-1- 


(0,0) 


(0,0) 


(0.380,0.250) 


mSP15 


1110 


760 


1097 


51.6 


-1- 


(0,0) 


(0,0) 


(0,0) 


mSP15 


1395 


554 


-175 


59.2 


+ 


(0,0) 


(-0.040,0.918) 


(0,0) 


mSP15 


905 


500 


1460 


54.8 


+ 


(0,0) 


(0,0) 


(-0.350,-0.260) 


mSP16 


520 


455 


620 


55.5 


+ 


(0,0) 


(0,0) 


(0,0) 


mSP16 


282 


464 


67 


43.2 


+ 


(0.912,-0.529) 


(0,0) 


(0,0) 


NUSP12 


2413 


454 


-2490 


48.0 


+ 


(0,0) 


(0,0) 


(-0.285,-0.848) 



Table 15.4: Benchmarks for the class HP where the Higgs boson [A, H) is the next 
nearest heavy particle after the LSP in mSUGRA and in NUSUGRA. The LSP and 
{A, H) sometimes are seen to switch. 
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Gluino Patterns (GPs) 



SUGRA 
Pattern 


mo 
(GeV) 


"ll/2 

(GeV) 


(GeV) 


tan /3 

(fu/Vd) 


M 
(sign) 


NUH 

(5h„, 5iij) 


NU3 

iSq3,StbR) 


NUG 

{5m2 ) ^Ms ) 


NUSP13 


2006 


1081 


-2027 


21.1 


-1- 


(0,0) 


(0,0) 


(0.207,-0.844) 


NUSP14 


3969 


1449 


-6806 


29.3 


+ 


(0,0) 


(0,0) 


(0.611,-0.834) 


NUSP15 


1387 


695 


2781 


50.5 


+ 


(0,0) 


(0,0) 


(0.136,-0.827) 



Table 15.5: Benchmarks for the class GP where the gluino g is the NLSP. Such a 
pattern was only seen to appear in NUSUGRA models with non universal gaugino 
masses. An analysis of light gluinos in the MSSM can be seen in [152] . 
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